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PREFACE 


This volume represents the second part of the final 
report for NASA Grant NSG-1414, Suppl. 9. The first part, 
designated Part I, was published in August 1987 and focused 
on the stability analysis of large space structure control 
systems with delayed input and the minimum time attitude 
slewing maneuver of a rigid spacecraft system with numerical 
examples based on the rigidized model of the Spacecraft Control 
Laboratory Experiment (SCOLE) orbiting configuration. 

This volume, designated as Part II, is based on the re- 
cently completed Ph.D. dissertation by Cheick Modibo Diarra, 
entitled, "On the Dynamics and Control of the Spacecraft Control 
Laboratory Experiment (SCOLE) Class of Offset Flexible Systems." 
First the open-loop dynamics of the orbiting SCOLE system are 
modeled to include the flexibility of the mast which connects 
the reflector to the Shuttle. The stability of system motion 
with respect to the nominal equilibrium during station keeping 
is considered for special cases, both for the 2-D and 3-D 
motion models. The control law synthesis is addressed for both 
small disturbances during station keeping operations and also 
during large amplitude preliminary slew maneuvers about the 
Shuttle's roll, pitch, and yaw axes, respectively. 
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Control law gains for both cases are based on linear quadratic 
regulator techniques. For the case of the rapid slew maneuvers 
the results presented here can form a basis of comparison with 
other results presented in Part I of this report and based on 
the application of two point boundary value problem techniques. 
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ABSTRACT 


A. mathematical model is developed to predict the 
dynamics of the proposed orbiting Spacecraft Control 
Laboratory Experiment during the station keeping phase. 
The Shuttle as well as the reflector are assumed to be 
rigid* the mast is flexible and is assumed to undergo 
elastic displacements very small as compared with its 
length. The equations of motion are derived using a 
Newton— Euler formulation. The model includes the effects 
of gravity, flexibility, and orbital dynamics. The 
control is assumed to be provided to the system through 
the Shuttle's three torquers, and through six actuators 
located by pairs at two points on the mast and at the mass 
center of the reflector. At each of the locations, an 
actuator acts parallel to the roll axis while the other 
one acts parallel to the pitch axis. The modal shape 
functions are derived using the fourth order beam 
equation. The generic mode equations ai a derived to 
account for the effects of the control forces on the modal 
shapes and frequencies. The equations are linearized 
about a nominal equilibrium position. When the interface 
point between- the mast and the reflector is assumed to 
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coincide with the mass center of the reflector, it is seen 
that the pitch equation is decoupled trum the roll and yaw 
motions. When the interface point is offset along the 
roll axis the pitch equation is still decoupled from the 
two other equations (roll and yaw). It is seen that the 
open loop system is unstable for both cases due to the 
(gravitationally) unfavorable moment .of ineertia 
distribution. When, in addition to the roll axis offset, 
a pitch axis offset is introduced into the system, the 
equations describing the roll, pitch, and yaw motions are 
seen to be all coupled together. It is further seen that, 
in the presence of gravity gradient torques in the system 
dynamics, the system assumes a new equilibrium position 
about which the equations will have to be linearized. The 
linear regulator theory is used to derive control laws for 
both the linear model of the rigidized SCOLE as well as 
that of the actual SCOLE including .the first four flexible 
modes. The control strategy previously derived for the 
linear model of the rigidized SCOLE is applied to the 

non-linear model of the same configuration of the system 

• * 

and preliminary single axis slewing maneuvers conducted. 

The results obtained confirm the applicability of the 
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intuitive and appealing two-stage control strategy which 
would slew the SCOLE system, as if rigid to its desired 
position and then concentrate on damping out the residual 
flexible motions. 
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CHAPTER ONE 


INTRODUCTION 

The problem of maneuvering a flexible spacecraft 
while suppressing the induced vibrations is becoming 
increasingly important. NASA is involved in studies which 
are concerned with the control of flexible bodies carried 
by a Shuttle in an Earth orbit. Similar experiments are 
being conducted in Earth-based laboratories. It is then 
desirable to derive a formulation which can accommodate 
both types of experiments. 

NASA is currently involved in at least two 
experimental programs to test techniques derived for 
active control of flexible space structures. 

In several versions of a recent paper, SCOLE^^ 
(Spacecraft Control Laboratory Experiment ) , Lawrence W. 
Taylor, Jr. and A.V. Balakrishnan have described the first 
which is ground based. It is a laboratory experiment 
based on a model of the Shuttle connected to a flexible 
beam with a reflecting grillage mounted at the end of the 
beam (Figure 1.1). As a part of the design challenge, 
the authors stressed **he need to directlv compare 
competing control design techniques and discussed the 
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feasibility of such a direct comparison. Concern would be 
given to modeling order reduction, fault management, 
stability, and dynamic systems. Ground-based 
experimentation has its limitations because it is almost 
impossible to to duplicate -the space environment in a 
laboratory. The second experimental program is known as 
Control of Flexible Spacecraft (COFS)^ 2 ^ and consists of 
experiments designed to control flexible bodies carried by 
a Shuttle in an Earth orbit. Because of the cost and 
risks involved in testing control techniques in space, 
COFS includes laboratory simulations of similar 
experiments which will precede the space test. Therefore, 
in assuring the success of both SCOLE and COFS, 
mathematical modeling and computer simulation are 
required. 

To accurately model and simulate flexible spacecraft, 

one needs a thorough knowledge of its structural behavior, 
f 3) 

In a paper v subsequent to the design challenge, the 
modal shapes and frequencies for the SCOLE system were 
derived. In the analysis of the mathematical model of 
reference 3, the SCOLE system is assumed to be described 
by partial differential equations in which the variables 
separate. The assumptions in that study did not create a 
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noticeable difference with the results previously 
derived^^. Based on the equations describing the motion 
of the SCOLE system provided in reference 1, the 
expression for the reflector line-of-sight (LOS) error was 
expanded analytically and studied c arefully . 

Analytical results shoved that the SCOLE* s LOS error is 
independent of the Euler yaw attitude angle so, only two, 
instead of, originally three, angular parameters were 
needed to be concerned with in designing the pointing slew 
maneuvers . 

Numerical simulation^ test results indicated, then, 
that the single axis bang— bang or bang— pause— bang slew 
maneuvers work fairly well for pointing the LOS of SCOLE. 
The best pointing accuracy and shortest slew time were 
attained when using the Shuttle torquers and actuators 
placed on the reflector while imposing a 5 degree/second 
slew rate limit on the design. 

Recently, a paper^^ concerned with the derivation of 
the equations of motion of the SCOLE class of flexible 
structures was published. The equations are supposed to 
describe a manuevering flexible spacecraft both in orbit 
and in an Earth based laboratory . The analysis is 
based on a perturbation technique in which the large 
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rigid-body motion are regarded as the unperturbed motion 
of the spacecraft while the induced elastic motions and • 
deviations from the nominal rigid motions are considered 
as perturbations. A maneuver force distribution *n the 
SCOLE system corresponding to the least amount of elastic 
deformation is derived. The paper axso highlights the 
coupling between the rigid and flexible modes. 

With the aforementioned papers as a background, the 
present study commenced by first reviewing literature 
pertaining to Reference 1, together with texts and papers 
which treat structural dynamics modeling and boundary 
value problems^\ 

Then, a mathematical model of the SCOLE system is 
developed assuming, that: the space Shuttle is a rigid 
body; the reflector mast is a flexible l)eam type 
appendage; and that the reflector is a rigid plate. The 
mast shape functions and frequencies are obtained from the 
fourth order flexural beam partial differential equation 
with different boundary conditions assumed to be imposed 
on both the Shuttle and reflector grillage ends. The 
system is represented as a beam connected at both ends to 
bodies with in-ertia. 
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Frequencies and modal shapes are derived for in-plane 
and out— of— plane bending modes as well as for the shaft 
torsional vibrational modes. 

The equation describing the in-plane dynamics of the 
system are developed based on an Eulerlan formulation. 

This equation is lineararized about a nominal motion where 
the Shuttle would have its velocity vector along the local 
horizontal. 

Also undertaken in this study is the modeling of the 
three dimensional dynamics of the SCOLE configuration 
based on the Eulerian techniques already employed, in the 
development of the in— plane (2-D) open-loop dynamics. The 
increased complexity of this three dimensional formulation 
should be emphasized. The techniques consist in isolating 
an elemental mass of the system in its deformed state and 

deriving it-s angular momentum taken at the mass center of 
the Orbiter. 

The position vector extending from the origin of the 
coordinate system to the elemental mass of the mast or the 
reflector accounts for the elastic displacements. The 
expressions for these displacements are derived from the 
mode shape functions generated during the three 
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dimensional structural analysis of the system (see 
Appendix A). 

The equations obtained for the elemental masses of 
the components of the system rn Integrated over the mass 
of the entire system to yield its angular momentum about 
the mass center of the Orbiter. The derivative of the 
system angular momentum vith respect to time is equated 
with the gravity gradient torques^ on the system about 
the same point (see Appendix C) . Such a vectorial 
equation, when projected along the three axes of rotation, 
yields the system rotational equations of motion. These 
rotational equations of motion are then linearized to 
yield a model which provides the basis for the control law 
synthesis developed in this study. But first, the 
stability analysis is conducted in the three following 
steps : 

first, when the mast is assumed to be rigid and to be 
connected to both end bodies, at their mass centers - ; 

second, when the interface point between the mast and 
the reflector is offset with respect to the reflector mass 
center in the "x" direction while the mast is still 
assumed rigid; and 
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finally when the "y" offset is introduced into the 

rigid SCOLE configuration. The flexible model of the 

SCOLE system, which Includes the effects of the first four 

flexible modes, is obtained by substituting the • 

expressions for the acceleration of the modal amplitudes, 

obtained from the generic mode equations, into the 

rotational equations of motion. These expressions take 

into account the effects of the controllers on the modal 

shapes and frequencies of the structure. 

( 9 ) 

The ORACLS v package is used to derive control laws 

for both the rigidized SCOLE (linearized model) and the 
linearized model of SCOLE with a flexible mast-. 

The control law, based on the linear regulator 
theory, derived for the linearized model of the rigidized 
SCOLE is also used for large amplitude rigid motion 
maneuvers. The non-linear equations describing the 
dynamics of that model are derived from the more general 
rotational equations of motion previously obtained. From 
the numerical results obtained from the simulations of all 
three models, conclusions are drawn regarding the modeling 
technique used herein, and the control efforts versus 
maneuver time of this strategy is compared with control 
laws previously presented for consideration for their 
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implementation in the laboratory model of the SCOLE. Under 

some additional assumptions, the equations describing the 

dynamics of the SCOLE system can be modified and adapted 

to systems^ with offsets proposed or currently under 

development such as the Wrap-Rib antenna^ 10 ^ in which the 

attachment of the lower mast, to the rest of the system, 

( 11 - 13 ) 

is offset and the tether connected platform (Kinetic 
Isolation Tether Experiment) in which the location of the 
interface point between the platform and the tether can 


vary. 
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DRAWING OF THE SCOLE CONFIGURATION 


FIGURE 1.1: 


CHAPTER TWO 


TWO DIMENSIONAL ANALYSIS OF THE SCOLE CONFIGURATION 


In this chapter, the equation describing the in-plane 
(rotational) dynamics of the SCOLE configuration (Figure 
II. 1) is derived using the Eulerian moment equation. 

The folllowing main assumptions are made in the 
development: 

a) the space Shuttle orbiter is assumed to be 
rigid; 

b) the mast, treated as a 130 ft long beam, is 
rigidly connected to the Shuttle at G, mass 
center of the orbiter, and to the reflector at 
0^, the interface point. 

c) the reflector is considered to be a flat and 
rigid plate with its mass center at G^; 

d) the mast is assumed homogeneous and isotropic; 

e) the elastic displacements are assumed small as 
compared with the length of the mast. 

f) the analysis here is performed by assuming that 
the mast vibrates at only one of its flexible 
modal frequencies but can be extended, within 
the linear range, by superimposing the effects 
of several modes. 

Angular M omentum of the Shnttle with Respect to its 
Center of Maas, G. 

The center of mass of the Shuttle is considered to 
move in an orbit about a fixed point, the geocenter 0; its 
angular momentum about its center of mass is given by 
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*S/G ' ( «-“o ) X S2 J " (II- 1) 

ft 

where 9 is the pitch rate, the angular velocity of the 
orbit, I g2 the moment of inertia of the orbiter with 
respect to an axis passing through 6 and perpendicular to 
the the orbit plane. 

II. 2 Angular Momentum, of the Mast With Reanect to G. the 
Orbiter Center of Mass. 

Consider here an element of the mast located at point 
P, with mass dm. The elemental angular momentum of such 
an element is given by 

dH M/G - GP x dJpP)^ dm (II. 2) 

dt 0 

where Rq is the inertial frame centered at the geocenter, 

0 . 

If one notes that 

GP - r Q + q; (II. 3) 

then. Equation (II. 2) may be expanded according to: 

Vg“ 

M 

il_ip is expressed using the relationship between the time 
dt |K 0 




(Fj+q) x (r^Kj)dn (H.4) 




12 


rate of a vector in an inertial (R Q ) and rotating ( R ± ) 
frames, i.e. 

4 I..--* I* +“*/* * * 

The in-plane transverse elastic displacement vector, *q~, 
(assuming a single mode of vibration) (see Appendix A) is 
found to be: 

T- cos (u»t+$) {A cos 8z+B sinBz+C coshBz+D sinhBz} i § (II. 6) 

After substitution of Equations (II. 5) and (II. 6) into 
Equation (II. 4) and integration term by term, one can 
develop: 

<‘ H M/a- ((r 0 +q)3t[fi S/R ^ 0 -K /EM +!l s/R0 >n[])dm 

_ + _ *► _ -*■ _ 

- < r o x( a S/R x V^/M^S/RO x ^^ x(n S/RO x? (P } ^ 

• * 

where one assumes q . q and q . q small as compared with 
the other terms. / This can be explicitly rewritten as: 

dfyg - { 6 +6 R -^z 4 ~z sinfofchfrMA cosBa+B sirSatC coshBaf D siriiBz} ]j dn 

with dm»6 dz where 5 is the mass per unit lengrh of the' 
mast, ^ S / R ^“ (® -w 0 )j» tlle inertial angular velocity of 
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the Shuttle , "5|j^the velocity of the mast element as seen 
from , the frame moving with the mast, and 


e * ’ at Ii ^ il «-i 




H M/G “ \ d h m/g 


H M/G “ ( e - w o ) p ~ + p sin (“*+♦)*) j 


where f - A (L sinBL + cosBL - 1 ) + B (L cosBL - sinBL) 

6 s 2 7 ~ ?r 


+ C (L sinhBL - coshBL + 1_ ) + D (sinhBL - L coshBL) 


2 2 
B B 


finally 

+ 2 

fyg" ( + M co ain^o t+4>) (f/L)]j (n.7) 


II « 3 * Angular Momentum of the Rigid Reflector About G. 
the Orbiter Maas Center "" 

The reflector being rigid, its angulai momentum can 

be found at G by a simple application of the transfer 

theorem (see Appendix B). 
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H 


R/G" H R|Gi +M R GG 1 x jd__ C GG i ) . R 

dt 1 0 


(II. 8) 


h r/g • ( ® + ®R’ (, ‘o )I R 2 J where ®R - i- [ la. ]| 

dt ' 2 «l 

0r— «8in ( wt+<f>) B[-A sinBL + B cosBL + C sinhBL + D coshBL] 
Because of the assumed magnitude of the transverse elastic 
displacement, it can be assumed that the length of the 
mast remains constant* Subsequently, 


dGG , | 
It -1 1 R. 


^n R + ywi* 


This can be expanded using Equation (II. 5) with the 
result 


4t 0G l!^ tL(e ’-“ 0 + ®R) 1 ‘n - 1 («-“o + V * M 
where X is the distance between 0^ and G^ 

GGj . -Lk M + Ii M 

G °l x (CG l>| R ■[(“•*■ \ -“o> L 2 + X 2 (0+e E -«) o )j] 

]j (II. 9) 


II. 4. 


Angular Momentum of the System About G. 


®syst/G " ± m i ®i/G (II. 10) 

® sjet/G m J S2 + («-“ 0 + ®R )ml2/3 

+ M u) sin (u t+4> ) f /L 

+ (8+VV (I R2 +M R [L 2 +I 2 ]))j (II. 11) 

RotationaJ^^gnation of Motion of the Open Loo p 
System • 

The rotational equation of motion is obtained by 
equating the rate of change of the angular momentum of thie 
system about G with the external torques acting on the 
system taken at G. Here it will be assumed that the only 
external torque is due to the gravity - gradient. (The 

effect of control torques will be treated later in Chapter 
Five) . 

f t ^syst/G 5 I * J - (11.12) 

which can be explicitly written as: 

e'a^+M^/S+I^+MjU 2 * l 2 ))+ e* R ci R2 + M l 2 /3 + Mr (X 2 + l 2 )) 
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2 

where -3 m o e (I^-I^) is the y component of the 
gravity-gradient torque acting on the system at its center 
of mass, S (see Appendix C). 


II. 6 


Linearization of the Rotational Equation* of 
Motion 

Let t , the dimensionless time, be equal to 


t. 


mi de de ' 

Then 37 "0 3T ' V 


and 


dt' 


2 d 2 e 2 „„ 

“0 7 2 * "0 8 ' 
a x 


Equation (11.13) can be rewritten as 


“ V " C3^2^ 2 /3 +I k -+M r (X 2 + l 2 )»\ 4f^(X 2 + L 2 )) 

-Mu 2 cos (o)t+*)f/L + 3w 2 (I 1 -I 3 )-0 (11.14) 

Dividing Equation (11.14) by ML 2 yields 
e"(I S 2 /ML 2 +1/3 + I R 2 /ML 2 +jp[(X/L ) 2 + 1]) 


+ 8 "r ( I R2 /Ml2 + 1/3 + m r t CX/L) 2 + 1]) 

M 

- (w /“ 0 ) 2 COS (o)t +<j>) f/L 3 + 3 9 (I 1 -I 3 )/ML 2 -0 (11.15) 


Equation (11.15) can be recast as 

c x e + c 2 e - f(t) (11.16) 

if u / Uq ■ W , M^/M “ w and X/L ■ X 


* - 
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then, C x - I S2 /ML 2 + 1/3 + I R2 /ML 2 + u (* 2 + D 
C 2 - 3(I X - I 3 )/ML 2 

f(x) - W 2 cos (wr+4> ) (f/L 3 ) - e" R (Ij^/ML 2 + 1/3 + (1 + X 2 )) 

11. 7. Stability Analysis of the In-plane Motion 

The stability analysis will be conducted in two 
phases: first in a torque free situation and second in 

the presence of gravity gradient torques. 

11. 7. A. Stability Analysis in the Torone-Frec 

Configuration 

In this case, Equation (11.16) becomes 

c x e " - f(x) 

which is integrated twice to yield 

0 (x) cos (Wx-te) (f A. 3 ) - e^i^/H 2 + 1/3 uCx 2 *!))*^-^] 

C 1 

The value of the constants Kq and are derived from the 
initial conditions (assuming 6 ■ 0) i.e. 
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when r . 0, 6(0) - e Q and e’(0) ■ 6^ 

C 1 e o " ( ’ f / L 3 ) “ 0 R (O) (I R2 /ML 2 +l/3+ u (a 2 + 1 ))+ 

or Z 1 - C x e 0 + e R (0) (I R2 /ML 2 +l/3+ (x 2 +1)) + (f/L 3 ) 

and 

c i e o " -^ R f C0)(i R2 7' MI ‘ 2 +i/3+ w(a 2 + 1 ))+K 0 

or K q - C^q' + e R ’(0) (I R2 /ML 2 +l/3 + (x 2 + 1 » 

finally, 6(t)- _1 { ( ( f / L 3 ) ( 1-cos(Wt ) ) + 0.(0. + 9 * ) 

c 1 10 0 

+ (I R2 /ML 2 +l/3+ y (X 2 +l)) ( 0 R (0 )t + e R (O)-0 R )} 
(11.17) 


A numerical simulation of Equation (11.17) (Figures 
11*2 - II. 6) for different values of the x offset and 
frequencies has shown that: 


a) the system oscillates about an equilibrium 
position, in the absence of gravity-gradient 
torques and disturbances, different from zero 
due to the forcing effects of the flexibility 
and the related coupling due to the offset. 
(Figure's II. 2 and II. 3) 

b) the amplitudes of the oscillations increase with 
the offset (Figure II. 4 and II. 6), according to 
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the coefficient of the second group of terms in 
Equation (11.17). 

c) finally, given an offset and an initial 

disturbance, the amplitude and the frequency of* 
oscillation about the equilibrium position 
increase with the frequency of the mode of 
vibration of the system (Figures II. 4 and II. 6) 
Note that the coefficients in f (A-D) are 
related to the modal amplitude functions. 


II. 7. B. Stability Analysis of the System In the Presence 

of Gravity Gradient Torque 

In the presence of the gravity-gradient torques, the 
two dimensional motion of the SCOLE system is described by 
Equation (11.16): 


c x 6 + C 2 « . f ( T ) 

In the absence of flexibility, f (t)*0, this equation 
reduces to 


Cj6 + c 2 e ■ 0 

Since C 2 ■ 3 (I^-I^/ML^ is negative for the SCOLE 
configuration, the solution, s(t), for this case is 
unstable. The case will be reconsidered in chapter IV 
where the three dimensional dynamics of the rigidized 
SCOLE is analyzed. 
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Here, Equation (11.16) is numerically integrated and 
the motion simulated for two different values of the 
offset (X- 18.75 ft and X- 37.5 ft): Figures II. 7, II. 8, 

and II. 9. 

In addition to the tendencies depicted earlier for 
the case of the torque free configuration, it is now seen 
that the SCOLE system, with the presence of 
gravity-gradient torques in its dynamics, is unstable. 
This is due to the inertia distribution of the system in 
the configuration considered here. 




SCOLE 2-D (NO GRAVITY-GRADIENT TORQUE) 






= SCOLE 2-0 (NO GRAVITY-GRADIENT TORQUE) 





SCOLE 2-D (NO GRAVITY-GRADIENT TORQUE) 






. SCOLE 2-0 (NO GRAVITY-GRADIENT TORQUE) 
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. SCOLE 2-D WITH QRAVITY-GRADIENT TORQUE) 





SCOLE 2-D (WITH GRAVITY-GRADIENT TORQUE) 
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CHAPTER THREii 


THREE DIMENSIONAL EQUATIONS OF MOTION : THE ACTUAL SCOLE 

CONFIGURATION 

In this chapter, the three dimensional formulation of 
the SCOLE dynamics is developed based on a Newton - 
Eulerian formulation. The Shuttle and the reflector are 
assumed to be rigid bodies and the mast is modelled as a 
connecting flexible beam. 

The expressions for the general displacement (See 
Appendix A) of an elemental mass on the mast are derived 
from the three dimensional mode shape functions consistent 
with the boundary conditions on the mass. The three 
dimensional rotational equations are obtained by taking 
the moment of all the external forces acting on each 
elemental mass, at some arbitrary point, and equating it 
with the moment, about the same point, of the inertial 
forces acting on the element. 

These equations must then be integrated over the 
entire system and then projected on the three axes of 
rotation in order to obtain the rotational equations of 
motion. Similarly, generic modal equations (See Appendix 
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D) for the flexible mast modes are obtained for the SCOLE 
system. 

III.l. SCOLE System Geometry 

Since the SCOLE system has three components with some 
relative degrees of freedom, it takes at least four 
coordinate systems to describe its geometry in its 
deformed state. 

Let therefore, Rq be an inertial frame centered at 

AAA 

the geocenter; (ij, kj ) , a frame connected and 

moving with the orbit with parallel to the angular 
momentum vector of the center of mass of the Shuttle and 

A AAA 

with ij directed along the positive orbit; R(i, j, k), a 

* 

frame centered at G, mass center of the orbiter; R2(i 2 » 

A A. 

$ 2 ' ^2^ a ^ rame ®oving with the reflector and centered at 
G, its center of mass. (See Figure III.l). 

If Mq is the orbital angular velocity of the Shuttle, 

then 

“o~Vi 


(in.i) 
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Let us assume the following Euler angle sequence 



then, ^ g/jj » the inertial angular velocity of the Orbiter 
can be expressed as 

S S/R„ - <*-V J 1 + * k ' + * i (III. 2) 

* t * A 

Since, k » cos k + sin <1* j and 

A A A A 

■ sin ♦ i + cos $ cos 'I' j - cos $ sin 4 > k, 

-► 

^S/R can be rewritten in the body frame as: 

fl s/R 0 - a x i + V + “ z k 

, A A 

■ [(e«-w q) sin $ +if» ]i + [ (0 — o)q) cos <j> cos if» '+ 4> sin i|/]j 

+ [<fr cos i/j - (0-Wq) cos <t> sin ij; Jk (III. 3) 

The reflector is assumed .to be rigidly connected to . 
the beam, its angular velocity is that of the end, 0^, of 
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the beam. The angular velocity, n R/ , s , of the reflector 
with respect to the Shuttle can be expressed as 



i + e Rjl 


+ *. R k 2 


(III. 4) 


where 


■ d_ 
R dt 


t-3v(z_,_t) ] 
3z 


| z--L 


0 - [0u(z_,_t)] | (III. 5) 

R dt 3z 1 2 L 


and i - d *(z, t) i T 
RdT R l z - L 

u(*#t) f v(z,t) , and 4> (z,t) are the in-plane, the 
out of plane, and the torsional bending mode shape 
functions of the beam, respectively. 

Assuming the following sequence in the beam motion 
relative to the.orbiter: 

i. Out of orbit plane bending; 

li. Bending in a plane parallel to the orbit plane; 

A 

iii. Torsion ‘about k 2 , 

the unit vectors in the intermediate coordinate systems 
are expressed as: 
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±2 - cos cos 6 R i + sin $ R cos ^ R j + sin $ R sin iji R k 

A ^ A 
J 2 " “ sin * R cos 0j i + (cos $ R cos tf» R -sin ^ sinSg sini|> R )j 

A 

+ (cos $ R sin + sin* R sin e R cos ip R )k 

* A 

kj ■ sin 0 R i - cos 0 R sin i|/ R j + cos e R cos i|^ k 

These relations can be recast in the following matrix 
format: 
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i 2 

a 

J 2 

" [T R*R2 J 

m « 

A 

i 

A 

j 

A 

k 2 

• 


T3T 1 

« 


(III. 6) 


in which the transformation matrix from the body frame, R, 
to the frame connected to the reflector, R 2 , has the 
following form: 


T Rh-R2T 


cos4^ cos^j sln^ cos^ 

«e<fe cos^ cos^-eLn^ sinU, siroj^ C<t^ St^fs±n4^ sin^ cosi^ 
sin^ -COS0JJ sirajig cos^ oos^ 
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Next, the rotational equations- of motion for the 
system will be derived by taking the time derivative of 
the angular momentum of the system at 6, the center of 
mass of the Shuttle, and by equating it to the external 
torques applied to the system. 

III«2. Angular Momentum of the SCOLE System 

III. A. Angular Momentum of the Shuttle About its 
Maas Center. G . ’ 

The angular momentum of the Shuttle, taken as a rigid 
body, about its center of mass, 6 is 

®S/G * h/G S S/80 ( m - 7 > 

(See Appendix C for the inertia tensors of the different 
components of the SCOLE system) . 

III.2.B. Angular Momentum of the Beam About G 

Consider an element of mass, dm, *of the beam located 
at some point, P, such that (Figure III. 1) 

GP - ** - T q + q (III. 8) 

where F Q ■ -zk is the position vector of P in the 

A A 

undeformed state; q (z,t) - u(z,t) i + v(z,t)j in which, u 
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and v are the x and y components of the mode shape vector 
(Appendix A). 


by : 


The angu lar momentum of dm about G, 


dH M/G is glve?1 


dH 


M/G 


r x 


d 

dt 


.(r) | 


R- 


dm 


where 


(HI. 9) 


. r » - 2 k + ui + vj 


(III. 10) 


Equation (III. 9) is expressed explicitly as: 

A A A A A A 

<®M/G ■ t (-zk+ui+vj)x d_ (-zk+ui+vj) | R ^} dm 

4LjR Q “dt^’^* r * ( u ”S 2 v -zfly)i+( v +0 z u+znx)j+(0x v “ u 0y) lc 

After substituting the different terms into Equation 
(HI* 10) the following expression results: 


dH M/G " {[*(^fl z u) + v(« x v-fi y u) — 2 2 fi x ]i 
+ [- z ( u— 0^ v ) + u(8 y u - « x v) + z 2 £2 y ] j 
+[u(v+ 0, u) - v(u-fl 2 v) + 2(ufi x +vfi y )]k}dm (III. 11) 
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where 

u(z,t) - l p n (t) S n (z) and v(z,t) - l p n (t) S° (z) 
n x n y 

Considering only a single node in the open loop situation, 
to show the form of various volume integrations, 
u - - “sin (“t+o) S x (z) and v - - a) sin (wt+Y) S y (z). 

Assuming small elastic displacements such that 

2 _ 

Si Sj/L <<1 and dividing by L , wh ere Wgis the 

Shuttle orbital angular velocity and L a reference length, 

then, 

^®M/G/ w 0^* - 1 ((zv+zft u+ fl z 2 )i + ( -zu 

2 z x 

“ 0 L 

+fl 2 zv+z 2 G y )j + (J 2 x uz + Q y zv)k}p dz 

where P is the mass per unit length of the beam. After 
“^tiplying both sides of this equation by l 2 , there 
results: 

dH M/G - {(*v+zu 8 -+z 2fl )i + (zu+zvfl +z 2 « )j 

“ z y 

A 

+(zu8 y +zv® )k) P dz 
x 7 


(III. 12) 
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The total angular momentum of the mast about G is 
obtained by integrating Equation (III. 12) over the total 
length of the mast. 


H» 


S 


-L 

dH 


l M/G " J an m/G (III. 13) 

0 

The eight terms appearing in dH M / G are integrated using 

integral tables - e.g. 

,-L -L 

dz- - -P <*>sin (wt+y) cosfiz+C^ sinhBzri^ coshBz) dz 

0 0 


J p w 


- P«sin (wt+y)[A 5 (^gin0L + L cosBL) + B, (L sinBL + cosBL - 1) 

2 7~ F 


+ C 0 (-L coshBL + sinhBL) + D„ (L sinhBL - coshBL - 1 ) 

6 e 2 6 "7“ 7 2 

To simplify the notations, let 


f^(B) " {A., L cos BL — sin BL ) + Bj (L sinBL + cosBL - 1) 

X Q O i n ~ 


2 2 
B . B 


+ C, (sinhBL - L coshBL) + D., (L sinhBL - coshBL + 1)1 

X o ft 1 . a ~ ■ _ 


6 2 6 2 


V 
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After substitution of f. and M for P, where M - mass 

1 L 

of the mast, in the expression of one arrives at: 

H m / g » M (P z cos (wt+a)f 1 - w sin («*> t+Y)f 2 -fi x L 3 /3}i 


+ [ a) sin (wt+a)f^+ & 2 cos (wt+Yjfj-fl L 3 ]j 


+[& x cos (a)t+o)fj+ Jly cos (wt+Y^Jk} (III. 14) 

III.2.C. Angular Momentum of the Reflector About. G 

Since small deflections are ^sumed for the beam, the 
reflector can be assumed to be located at a constant 
distance from 6, the Shuttle mass center. 

Using the transfer theorem for the angular momentum 
(See Appendix B), the angular momentum, H r ^q, of the 
reflector, assumed rigid, about G can be expressed as: 

4/G * ^ n R/E 0 + Wr <*1 * ij. (HI-15) 


where the inertia tensor of the reflector expressed 


1 R/Rq “ fi R/S +fl S / R 0 


at Gj , its center of mass, and fl, 

(respectively, the inertial angular velocit-y of the 
reflector , its relative angular velocity with respect to 
the Shuttle and the inertial angular velocity of the 
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Orbiter) are both expressed in the same coordinate system 

* moving with the reflector. 

The inertial angular velocity of the reflector is 

expressed in R 2 using the transformation matrix [Tij] DH . D 
^ K K 

as : 


t 


2 


a R/E 0 - {a j + i„) T n + a y T J2 + a 2 t 13 + s R sin » R )i 2 

+ ((n x + *„> T 21 + a y t 22 + a z t 2J + e R cos * R )j 2 

+ << n x + *R> T 31 + a , T 32 + T 33 + V*2 

A , A A 

" Q 1 i 2 + a 2 J 2 + fl 3 k 2 (III. 16) 

Now, after rewriting the second term in Equation (III. 15), 

Mr <*! * i- - MgCGO.+O.G.) x d_ (GO,+o7g.)i r 

dt R. dt A 1 


where 0^ is the reflector attachment point to the mast, 
Tt lg " + v(-L,t)J) + SJgyjj x (-Uof uL + vj)!^ 


-(u-G y L - 
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JT°l G l |* 0 “°R|R x °1 G 1 - <n l i 2 +a 2 J 2 +a 3 it 2 ) 1 (X h + Y ^ 

* (fl j ^ 2 "* 83X j 2 — ftgY ±2 
1, and Y are the ”x w and ” 7 " offset coordinates, 
respectively. 

After subsitution of the terms into Equation 
(III. 15), one arrives at: 

■f ^ 

H R/G *°1 I R 1 *2 + n 2 X R2 ^2 + *R3 k 2 

+ Mg { (bL + c(v+Y))i-(aL+c(u+X))j + (b(u+X)-a(v+Y))k) (III. 17) 

^Rl» Ijj 2 » an ^ !r 3 are the principal moments of inertia of 
the deflector. 

a - (u- fl y L -a £ v- « 3 Y T 11 +X0 3 T 21 +CQ 1 Y-0 2 X)T 31 } 
b - (v+n x L +n z u-n 3 Y t 12 + ft 3 x T 22 +(n 1 Y-« 2 X)T 32 } 

and 

c - <°* T -° y « - °3 1 T 13 + n 3 1 T 23 + (t! l T * V )T 33> 
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III. 2. 9. Angular Momentum of the System About G 


The angular momentum of the system about G 
is given by the sum of the angular momentum of 
three components evaluated about the same point 


* H syst/G ’ 
each of the 

* G . 


H syst/G " H S/G + H M/G + H R/G 

The unit vectors i 2 » j 2 » and k 2 are transformed into unit 
vectors along the Shuttle axes as follows: 

A A A ^ 

i 2 -cos* R cos0 R i + sin* R cos^ R j + sin* R sint|/ R 1c 

h " ~ shl \ coee R 1 + cos ^ R - sira^) j + (sin^ sinegcos^ 

A 

+ costjj sini|^)k 

^2 ™ 3in0 R i — cos 6 R sin^ R j + cos0 R cos* R k 

After this substitution into Equation (III. 17), the 
angular momentum of the system is expressed as: 

^yst/G " +£l as (<-»t+a)f^ - sin (u>t+r)f 2 

3 

+ & x L/3 ] (bL-tc(wy» +n x cos * R cos 9 R -Q 2 3^ sin* R cose R 

A 

+ n 3 ] K aii.e R ) 1 + 6 Lj + m [a) sin (cot+oOfj 

L 



\ 
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+ n 2 008 ] - ^ (aL 4c(ufX))+ Qj Ij^ sin 4> R cos4^ 

+ a 2 Igr (coe * R oo^ - sin^ sin e R sia^) -^3 Ijg cos ^ sin^Jj 

+ { - Q t 3 S 4 +Q z ^33 +^[fl x cos (ojWa)fj + fly cosfot+Y^] 

+ (b(ufX) - a(vfY» + Gj 3 ^ sin $ R sin'i'p 

+ n 2 ho. (sin 1? COB \ + °» $ R s^V + n 3 Ijq cos6 r oo^ R }k (m.18) 
or 

* V + V + V (III - 19) 


m*3« Rotational Equations of Motion (Torque free) , 

The rotational equations of motion for the system, 

when free of all external torques, are obtained as: 

• . 

^ H syst/G^R 0 . * H syst/G/S + n S/R 0 x H syst/G " 0 ( III * 2 °) 

The vector equation (III. 20) itself is equivalent to 


H + g H 

x y z 

H + 0 H 

7 2 * 

H + ^ H 

z x y 


0 

z 

0 

x 



0 H > 
y * 


(III. 21) 


Under the small angle approximation assumption, on 
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, 0, and 4> (sin a^a ; cos cr^l) the torque free rotational 

equations of motion (III. 21) for the SCOLE system, after 
linearization are rewritten as: 

i) The Roll Eauation: H + fl H - n n - n 

x y z z y 

(i-WQ^^si - (++V>I S4 + M [(♦+ * 0 h cos (o>t+o)f 1 
— u»( i^f ( i)q(|> ) sin (wt+a) fj - w 2 cos (cot-Hy)^ + (i-w 0 ^)L 3 /3] 

+ *+ ir®b*R )] ®l + “t)*R J R2 • (w 0* ‘V }I S4 

• 2 9 

+(u Q i+u) 0 t|i) ( I S2“ I S3 +I R2~ I R3 I Rp + - K“o*-“d^ cos ( wt+a ) f i 

3 

+(u, 0 *Ho 2 <|/) 3-+ ( W 2 -2(^0) cos (a»t+T) f 2 ]^ R I R2 

“<V_ 0 + “oV^ + “o^V^Ri) + Mr( u + X > — 0 )v 

+IXu>q<|>-w 0 $) - (a) 0 ^+WQt|»)u - X (w 2 iJJr + u> 0 * r ) 

2 

+ a^) i|»R X] - M^v + Y)[ (6- w () )u + ((O o 0-U) 2 )L 

+ ( w o^, w o*) v + Y (“0% + “oVV +w o*> - u o 0 R X 1 
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•• ••• B(i 

+J^L[v + + X + $ R $ R )] + Mg vX U g 

•• •• 

4^ Y(Y(H tf> R -^-a) 0 ^) - X (0+65)} - 0 (III. 22) 

ii) The Pitch Equation : H y + n z -jj x - 0 

•• 2 •• • 

0I S2 + ^ f“ “ cos (wt+a)^ + (<p+<yli) cos (wt+Y)f 2 

-2^C++<o 0 4t) sin (wt+r) f 2 + 6L 3 /3 + U(J <$ -w Q <t>) cos (wt+r)(f 2 )] 

+ (®*^) Ir 2 “ [u “(&+6 r )L- Y ($+$ R + u^ + wq^)+ 

6r X] - 2 Mjj Xu^q u - Mr XY Ql>- a^+ij; R -Wq$ r ) 

+ Mr X 2 (e+€^) - Mr Y (*-a> 0 <j>) ojqL- 0 (III. 23) 

iii) The Tew Rogation s H z ®y ” " 0 

“(^ fu to *S4 + (♦ +w o^ ^S3 f X R3^ + — [(4>“wg$) cos (“ ,t+a ) 

L 

-0)( sin (Mwa)^ +Qcae (u>W)f 2 -u(^ sin (u>t+Y)f 2 ] + ^(irtX) fr-J/L -a» 0 *L 

+ ♦ u + X ( “o *R + * R + “0 f )i + * u t^ <$“ < ^0 *- ) L + ( ♦ + W 0^ ) U 



+ X (4 Uq ’H«^+w 0 <J' r )] ” ® L “ * v “♦ *“*“ 0 *" ojq ip v-ipu> Q ▼ - Y ( V R + 

+ * r +“ 0 4 ♦) + Xe R ] - ^ v[u - (e- o)q)L 

-<4«fo&v +u 0 e R x - T (♦ +“ 0 <* w Ql 

3 

- ( 4 u k < fr) ,l ^j 1^2 + — ( ♦— u>q$) (u sh* ( wt+a )f^ + oj^ ) 

L 3 

-U) 0 (i-tOQt) (*-^ +> Cu+u»(f) 

-J^(ttriD^(* -ajQ^) + “o *S1 +w 0 ^ +u 0^ *S4 

+ M [o b (*+ Wo i(»)co8Ca)t+a)f 1 + (e- u, 0 ) “sin ( «t +Y)y 

L u 

+ utf 4a» 0 <H- ^r« 0 * R ) % -Hoq *r !r2 

i. + 0+“ Q +) L 2 + XL (♦ R U) Q + * r + uq + i) -t^XL ij, R 

+^Yu + Y^C- ♦ + ♦ jj^dyniQ^) + ? X (0f 0 -p} ■ 0 (111.24) 

Equations (III. 22 - III. H4) have been developed 
under the assumption that only a single mast flexible mod 


is excited* 
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For multiple mode interaction, terms involving 
cos (ut+a), sin (wt + a), etc. would be expanded to 
Include the effects of the multiple frequencies. Where u, 
v and their derivatives appear explicitly, multiple modes 
can be included by direct substitution. 





CHAPTER POUR 


STABILITY OF THE SCOLE SYSTEM IN SOME OP ITS 
CONFIGURATIONS 

In Chapter III, Equations (III. 22), (III. 23), and 
(III. 24) describe the dynamics of the orbiting SCOLE 
configuration. In what follows, the stability analysis of 
the SCOLE system will be conducted in three different 
steps. 

First, it will be assumed that the mast is rigid; 
also that the interface point between the beam and the 
reflector is the reflector center of mass; second, still 
assuming, the mast rigid, the interface point will be 

In the "x" direction; finally, a two dimensional 
offset of the interface point will be introduced. The 
mast will still be assumed rigid. The system dynamics, in 
all the aforementioned cases, includes the gravity - 
gradient torques. (Appendix C) 

IV - The SCOLE System Without Offset or FexibilitT 

In the absence of flexibility (f.-f 0 - iji « \p - i|>» 0 - 

i 4 K R R R 

t •• •* 

®r " h m V V V» tt *u-u-v-v-v-0) . and in the absence of 
offset in the location of the interface point (X-Y-Q),. 
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equations (III. 22), (III. 23), and (III. 24) respectively, 
can be rewritten in the presence of gravity gradient 
torques as: (See Appendix C) 

i|>[I sl + ML /3 + M r L 2 + I R1 ] - <fr’l g4 - o, 0 i[I sl - I g2 + I g3 

+ ^Kl “ *32 + W “ w 0**84 “ “ 0^ ^*S3 “ *S2 + *R3 “ *R2 ” 

-M r L 2 + 3(I 3 - I 2 )] - 0 (IV. 1) 

e[Is2 + + + + 3 <d 2 9 (I^ - I 3 ) + 3oiqI 4 «0 (IV. 2) 

-il S 4 +i(I S3 + I R3 ) + V [I S1 +I S3’ I S2 +I R1 +I R3” I R2 ]+ 

-“0 [I Sr I S2 +I Rr I R2 ] ^o^ X S 4 + 3 I 4 ] - 0 < IV * 3 > 

It is seen that in such configuration, in the linear 
range, the equation describing the pitch motion (Equation 
IV. 2) of the system decouples from the equations 
describing the motion in the two remaining degrees of 
freedom (Equations. (IV.l)and (IV. 3)). 
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IV , I • A • Stability Analysis 

Equation (IV. 2) can be recast in the following form: 
•• 

ehj - 0h 2 + h 3 - 0 (IV. 4) 

in which. 

h l - X S2 + X R2 + M R l2 + iSt 2 

h 2 " ~ I I^ and h 3* " 3w o X 4 

The homogeneous part of Equation (IV. 4) yields the 
following solution: 

6t -fit 

0 h - C x e + C 2 e (IV. 5) 


where s *\f h 2 /h. - 0.00176 based on 

nominal SC(j)Lg 
parameters * 

since for this configuration, h 2 /h 1 >.0, e h (t) is 
unstable. A particular solution to (IV. 4) is obtained as: 

0 • h-/h ? » 0.0012 (based on nominal SC0LE 

p .. parameters) (IV. 6) 

The constants of integration, and C 2 , are determined 

from the initial conditions as: 
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and 


C 


2 




giving 

0(t) ■ 6j+0p« 0 q cosh~ +^adrfiS t + hg/hjCl-oosh 5 1) (IV. 7) 


In the absence of control, the system is seen to be 
unstable in its pitch degree of freedom. 

Equations (IV. 1) and (IV. 3) which have the following 
forms, respectively. 


• • •• 

♦ kj + * k 2 - ♦‘kg + 6 k 4 - *k 5 - 0 (IV. 8) 


where; k^ 


k 


5 


+ ^ng - $n^ + « 0 (IV. 9) 

- [I S1 + ML 2 /3 + M r L 2 + I R1 ] { k 2 - -I S4 ; 

2 

“0^ I Sl“ I S2 +I S3 +I Rl” I R2 +I R3^ ’ k 4"“ w 0 I S4 ; 
“0 [I S3- I S2 +1 R3- I R2-( Ml2/3 > +3(I 3 -I 2 )-M i! L 2 ] ! 


n l*“ I S4 ; n 2"* I S3 +I R3 ; n 3 =a) °- I S1 +I S3 _I S2 +I R1 +I R3” I R2 ^ 
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n 4-“o tI sr I S2 +I Rr I R2 1 ! “ d "5-"o (I S4 +3I 4 ) 

can be recast in the following state matrix format: 


X ■ A X or 
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Some of the eigenvalues of the state matrix, in this 
subcase have positive real parts, based on the actual 
SCOLE system parameters indicating instability in 'the open 
loop dynamics of the roll and yaw degrees of freedom. 

IV. 2. The SCOLB System Without Flexibility but With 
Offset in the "I" Direction 

The configuration analyzed in section (IV. 1) is 

upgraded to the one considered here by letting X be 

non-zero in the equations of motion (III. 22) - (III. 24) 

and by setting all the flexibility terms and the "Y 

. offset" equal to zero. 

The equations of motion then become: 

i[I sl+ ( ML 2 /3) ^L 2 + I r1 ] - Jcig^XL) + I S3 - I S2 

+I R1 +I R3 _I R2 ] " “o* (I S 4 +M R IL) o't , [I S3 -Is2-iSi 2 +I |t3 

” I R2" M R l 2+3(I 3" I 2 vJ " 0 (IV. 11) 

• » 

9 IS? 3 - l') + 0 (1V.12) 

-tflg^XL] +i[I S3+ I R3 +M R X 2 ] ♦ Vt I sr I S2 +I S3 +I Rl 
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+I R3” I R2 “o' ^ I Sl~ I S 2 +I Rl“ I R2“ t % x2 ^ 

- u» 2 * {3(I S4 +M r XL)J - 0 (IV. 13) 

Again, it is seen that in this configuration, the 
pitch equation, (Equation IV. 12), decouples from the roll, 
(Equation (IV. II)) and yaw (Equation (IV. 13)) equations 
and can be rewritten as: 

eh^ -ehj* + h^ - 0 (IV. 14) 

where, h^ - *s2 +I R2 +M R (x2+ L 2 )+ML 2 /3 

b 2 - 3u 0 2 (I 3 ' - 1^); and hJ-Srn^ I g4 +M R X L) 

t f 

Here again, hj/h^ is a positive -quantity. By analogy with 
the previous configuration, 

e(t) - cosh 6 *t +-^sinh 6t + h^ /h '( 1-cosh St) (IV. 15) 

with s m Vh /h m 0.00176 (based on SCOLE. 

nominal parameters *°) 

In the absence of control, it is seen that the pitch angle 
• • 

is unbounded indicating an instability in that degree of 
freedom. 
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A reasoning similar to the one previously done for 
the case without offset, enables one to recast Equations 
(IV. 11) and (IV. 13) in the following state matrix format: 


where, p 


X 

- 

A X 

or 



» 






* 


0 

0 

1 0 


* 

• 

♦ 

. 

0 

0 

0 1 


♦ 



I 

1 

1 1 



i 


P 2 

P 4 

P 1 P 3 



•• 


1 

1 

1 1 


• 

♦ 


- p 6 

P 8 

-P 5 -P 7 \ 


♦ 


1 

1 


1 1 1 

I 

i 

1 

k 0 

n 0 


, n,k*+k_. 

I- 


2 

3 

■m • n 

' 5 2 5 

2 . 


II Ilf 

k l n 2" n l k 2 


• t 

k . n. 
3 2 


’ 2 ~ 1 * * » » 
k l n 2 _n l k 2 


(IV. 16) 


II II 

k l n 2-"l k 2 


II If 

-(n 4 k 2 -k 4 n 2 

I * ft 

k l n 2~ n l k 2 


Iff III 

i 2 +n 3 * 

p 5 - — ; p 6 - — 

n 2 n 2 


I 



n l P 3 
? 


• ~ n l p 4+ n 4 
p 8 


I 



- ^ 2 +m r l 2 +i ri' > 

■ “0^ I Sl +I S3“ I S2 +I Rl +I R3' I R2^ i k 4* "o ^ I S4 +M R n, ^ : 
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k 5-“o [I S3' I S2 +I R3- I R2- M R L2 -l^ 2 +3(1 V^], 
n l “ ’ n 2” I S3 +I R3 +M I! j2 ’ 

t 

n 3" “0^ I S1 +I S3" I S2 +I R1 +I R3” I R2^ 
n 4* a> o^ I Sl” I S2 +I Rl” I R2”^R^ 2 ^ an< * 

* 2 

“ -«Q{3(Ig^+MjXL) } . Here again, some of the 
eigenvalues of the state matrix have positive real parts. 
Therefore , the open loop dynamics of the system are seen 
to be unstable in its roll and yaw degrees of freedom. 

IV * 3- T he SCOLE System with Offset in Both the "I" and 
"Y" Directions but Without FlexibilitT 

If once more the description of the system dynamics 

is upgraded by introducing the "Y offset”, the rotational 

equations of motion become: 

* [I S1 +I R1* y * (l2+t2)1 *♦ X xz * *“0 [I S1 +I S3 

+2M R y2 - I S2 +1 Rl +I R3- I R 2 ) * “oS 11 -“on i S3- I S2 +I R3- I R2 
+ M R (Y 2 -L 2 )-W ; 2 + 3(I 2z -I yy )]- u 2 0 *I XJ -3 l , 2 o9 I ]ty 
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+ w ^ m r Y L + 3I yz ] - 0 (IV. 17) 

e[I S2 +ML 2 +Ir 2 + m r(L 2 +X 2 )] + *MrYL + ip M r XY+ ojq $ M R XY 
-w^MrYL+o^MrYL - 3u, 2 *I x "+3 u 2 9(I XX -I„) 

+3u, 0 I xz " 0 (IV. 18) 

♦ CI S3 +I R3 +M R (x2+i2)] - ♦ I zx +u b' ,; t I Sl +I S3" I S2 +I Rl +I R3“ I R2 
+2 M r Y 2 ] -0 MrYL - « 0 0MRXY + 3 W 2ei yz - w 2 K 4I x2 ) 

-u, 2 ^ i Si- i s 2 + i ri- i r2+ m r(Y 2 -X 2 )] - u> 2 M r XY-0 (IV.19) 

It should be noted here that the pitch equation no 
longer decouples froo the roll and yaw equations. 
Equations (IV. 17), (IV. 18), and (IV.19) can be recast in 
the following state matrix format 
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. »» 

X « A X + C or 


♦ 


0 

0 

• 

0 

1 

0 

0 


' 

<P 


0 

• 

0 


0 

0 

0 

0 

1 

0 


6 


0 

• 

* 

- 

0 

0 

0 

0 

0 

1 


♦ 

+ 

0 

•« 

• • 


a l 

a 2 

a 3 

a 4 

a 5 

a 6 


* 


a 19 

8 


a 7 

a 8 

a 9 

a 10 

a ll 

a 12 


• 

e 


a 20 

• • 

. ♦. 


. a I 3 

a 14 

a 15 

a 16 

a 17 

a 18 


• 

<fr 

. 

1 

a 21 


(IV. 20 ) 


where, given 

I xx- i si +i hi^ 2 +m r( l2+1 ' 2 ) ‘ Wy 

I yy -I S2 +I R2^ 2+M R( 1 ' 2 + ^ + M R JtL 

I zz- I S3 +I R3 +M R( x2+r2 >* V*"* 11, 

* . M_XY , 

A l -V 1 **- *2-T~ ■■ 

XX 

K i ■“0 M R rL/ ’ 1 n : *5 -4 “ 2;i 22- I 7 y )/I xi : 

V“o Z xz n xx ! A 7* 3 “0 1 xj /1zx> A 8“ , 0 (M R TL+3I jr 2 )/I x* 
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Here A^-A^t B^-B^ have been identified with primes in 
order to avoid confusion with coefficients appearing in 
the spatially dependent functions s (z), s (z) and e (z) 

* y 

appearing in the solutions to the partial differential 
equations describing the vibrations of the mast, Appendix 
A) . 


m r yl m r xy m 

A 9 " T~ ; A io“ I ; A 11 ‘ “oi 


MrXY 

w n i — ; A io “ 


^ ! A 14 ■ 3 “0 I xy /I yy* A 15 - ’“o'V 1 - 1 ' 1 . 


13 0 I 


zz yy ; 


M YL 

Al6"3“ 0 I xz^ I yy* A 17 " I zx / ^ I zz’ A 18 " I * 

zz 


M XY 

A 19" u 0 ^xx + *zz ” I yy^ I zz ’ A 20 “ w 0 I ’ 

zz 


A 21' 3 “ 2 V 1 .. ! A 22- 4 “ 2 l xz n zz 


A 23-"^ I xx- I yy- + ^z^hz 1 A 24-“X XY/I zx 


♦ f » 


®l"^ A 2 +A l A 18^^“ A 1 A 17 ^ j B 2*” A 1 A l 9^ ^ 1 ~ A 1 A 17^ ; 
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B 3* A 3^^ 1 “ A 1 A 17^ ; B 4*^ A 4 _A I A 20^ /(1_A 1 A 17^ ; 

B 5 -( a 5 + a i A 22)/(1- a i a 17 ); B 6 - (A g +Aj A 23 )/(1-Aj A j7 ); 
B 7 «( a 7 - a i A 21 )/(l-A 1 A 17 );B g - (Aj A 24 -Ag) / ( 1-A| Aj 7 ); 

Bg-(1+A 9 A 18 )/(A 10 +A 9 A 17 ); b io*^ a 12 +A 9 A 19^ A 10 +A 9 

a i 7 ); 

B ll“” A ll^ A 10 +A 9 A 17 ^ j B 12* A 9 A 20^^ A 10 +A 9 A 17^ ; 

B 13*^ A 14“ A 9 A 22 )/(A 10 +A 9 A 17 ); B 14"“ (A 13 +A 9 A 23 /(A 10 +A 9 

a 1 7 ) 

B 15»( a 9 A 21 -A 15 )/(A 10 +A g A 1? ); b 16 — C a 16 + a 9 A 24 )/(A 10 +A 9 

a 17 ); 

B - B 

*20- 4* " 9 a * *1 - *13 + *9 *7’ *2 * *9 *8 + B 15 ! 

* B 


a 3* B 9 a 9 + B 14 ’ a i 


B 9 a 20 + B 10 ; a l 


B 9 a ll + B 12 ’ 
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a 6“ B 9 a 12 + ®11 ! a 19 " B 9 a 20 + B 16 ! 


a 7 - ' B 5 

B 1 ’ B 9 


a 


- B 15 ~ B 7 

of > 

8 B l- B , . 


a 


. ~ B 6 
9 “ 1 

B 2 - B 9 


a 


10 


B 10 + B 2 


'ii 


B 12 * B 4 

b ; - B, 


12 


B 11 ~ B 3 
B i - B 9 


l 13 " A 17 a l + A 18 a 7 + A 22 ’ a 14 * A 17 a 2 + A 18 a 8 “ A 21 ; 


a 15 " A 17 a 3 + A 18 a 9 “ A 23 ’ A 16 " A 17 a 4 + A 18 a 10 “ A 19 


a 17 " A 17 a 5 + A 18 a ll “ A 20 * a 18 " A 17 a 6 + A 18 a 12 ; 

and a 2i“( A i7 a i 9 + A ig a 20 + A 24^ 

Since the Shuttle axes do not correspond to the 
principal axes of the system, the system dynamics appear 
in the following state form: 



indicating that the system equilibrium position ‘is no 
longer - ®0 " *0 " °* 

Let \ p , e . and a be the equilibrium position for 
this configuration of the system. Then, 

♦ " <l» e + n l and ♦ " ni 

0 • e e +n 2 and 9 " n 2 

♦ ■ *e + *3 and n 3 

The new state vector is [ , n 2 » 

Also * , 0 , and a satisfy 

c c 8 


a l 

♦e 

+ 

a 2 

®e 

+ 

a 3 

^ e 

■“•19 

a 7 

^e 

+ 

a 8 0 

®e 

+ 

a 9 

♦e 

" _a 20 

a 13 

*e 

+ 

a 14 

®e 

+ 

. a 15 

♦e 

«- a 2i 


• • • iT 

n 3 * nj * T12 • J * 


this simultaneous system is solved using 
[a] » [A] [ 6 e , ♦ e ]a-4- [^ e » e ft » * 


[A' 1 ] [a] 
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* - - . 548552 rad. 
e 

® ■ - . 019345 rad. based on nominal SCOLE parameters 

♦ ■ . 207672 rad. 

B ■ — 


For the open-loop dynamics of the non-linear model of 

the SCOLE system, the envelopes of the Euler angles are 

. »_ 

depicted in Figure IV. 1, as a function of time. 

After substituting the new state vector in the 


equations describing the system dynamics, linearizing them 
about the new equilibrium position, recasting them into a 
state format, one arrives at a system which can be cast in 


the following form 


(where a a ^ g are 


constants) 



(IV. 22) 


The open loop system in this configuration is also 
unstable due to the unfavorable intertia distribution. 







CHAPTER FIVE 


CONTROL SYNTHESIS 

In this chapter, the different components for the 
two-stage c ontrol strategy which would slew the system as 
if rigid and then damp out the mast vibration will be 
analyzed in the following manner: 

first, within the linear range, the motion of the 
rigidized SCOLE is controlled using a strategy, based on 
the linear regulator problem when the system is subjected 
to some small perturbations in its degrees of freedom; 

second, still within the linear range, the motion of 
the actual SCOLE system, including its first four 
vibrational modes, is controlled using a control law based 
again on the linear regulator theory when the system is 
subjected to initial perturbations in its different 
degrees of freedom. 

third, the control strategy derived for the linear 
model of the rigidized SCOLE is applied to the non-linear 
model of the same configuration. Preliminary slew 
maneuvers are tested by assuming single axis initial 
perturbations of 20° in the roll, pitch, and yaw degrees 
freedom, respectively. The three Shuttle torquers and 
the two actuators on the reflector are then assumed to be 
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the only sources of control moments. The controllers are 
seen not to reach saturation. 

?.l. Control of the Rlnldined SCOLB 

During the control of this model. It Is assumed that 
the actuators located on the mast (proof masses) are not 
activated. As a result, the system is controlled by means 
of the Orbiter torquers and the actuators located on the 
reflector (Figure 1.1). 

Since the Shuttle is equipped with three torquers 
acting about the x, y, and z directions, the total control 
torque available can be written as 

T " {M x U x + 130F v y> * + (M U y “ 130 F x v x )i 

+ < M 2 U z + 32 - 5 F x v x + 18.75 F y v ) k) ft.-lb (V.l) 

with the limits for M , M and M - 10,000. ft. lb: F 

* y z x 

and F y * 800 lb.^ The constraints, therefore, are 
l u x li. I u y l <. ll l u 2 l <. 1» l v x l < li and |v y | < 1 

where U, the control vector is expressed as 
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U ■ [v^» Vj, U lt t^. 0 3 ] , while the control 
influence matrix can then be written as: 


B 


0 

0 

0 

0 

-130F 3 
32. 5F 


0 

0 

0 

130F 

7 

0 

18.75F 


0 

0 

0 

M , 

-0 

0 


0 

0 

0 

0 

M 


0 

0 

0 

0 

0 

M 


(V. 2) 


The optimal control U which minimizes a performance 


index 


J ■ 


S 


(X T QX+U T RU) dt 


is given by 


U - -KX - -(R“ 1 B T P)X 


(V.3) 


where P is the positive definite solution of the steady 

Q 

state Ricatti matrix equation . 

The equations describing the closed loop system can be 
recast in the following matrix format: 


X - AX + BV 
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After substitution of -XX for U, the closed loop equation 
can be rewritten as 

e 

X - (A-BK)X ( V . 4) 

A parametric study was conducted by first examining 
the variation of the real part of the least damped mode as 
a function of different values for the (assumed) diagonal 
Q and R weighting elements (Figures V.l and V.2). In 
this initial study, each of the diagonal Q elements were 
assumed equal i.e. Q»diag. [SQ] and also each of the 
diagonal S elements were assumed equal R ■ diag. [SR]. 
Figure V. 1 . corresponds to a model of the rigidized SCOLE 
system where the dimensionality of the state vector is 6 
x 1 and 3 Shuttle torquers plus 2 reflector actuators 
describe the control inputs. On the other hand, Figure 
V.2 corresponds to the model of SCOLE including the first 
four flexible mast modes. For this case the state vector 
has dimensionality 14 x 1 , and four additional control 
actuators are assumed to be placed on the mast - two 8 ? 

1/3 the total length and the remaining two at 2/3 the 
total length (see Figure 1.1). 
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It can be seen from both Figures V.l and V.2 that the 
best closed-loop transient results are obtained from using 
larger values of the state penalty along with smaller 
values of the control penalty elements. However, when 
the closed loop dynamic responses were simulated using the 
best combinations -of Q and R it was seen that some of the 
controllers reached saturation levels for responses with 
initial conditions on pitch, roll, and yaw taken within 
the slewing angle range (i.e. approx. 0.3 rad.). 

As an alternative, the concept of split weighting of 
both the state and control penalty elements was 
considered, initially for the rigidized SCOLE model. 

Since the roll (and to some extent also the pitch) are 
easier motions to excite than the yaw, due to the SCOLE 
moment of inertia distribution, it seems intuitively 
correct to relax the penalty of these control inputs as 
contrasted with the remaining control penalty elements. 
Also since both position and rate feedback of the Shuttle 
rotational motion will be utilized, it appears logical to 
place a far greater penalty on the (angular) position 
displacements. Based on this philosophy and by trial and 
error t the set of Q and R which produced the largest 
absolute value of the real part of the least damped mode 
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(while at the same time avoiding saturation during 20° 
single axis slewing maneuvers) was selected as 

Q - diag. [5xl0 12 , 5xl0 12 , 5xl0 12 , 1, 1, 1] 

and R - diag.[l, 1, .1, .2, 1] 

For this set of Q and R thee closed log eigenvalues for 
the rigidized SCOLE model are calculated 
RCXj) Im(X i ) 

-0.431436E+02 
-0.431436E+02 
-0. 132023E+03 
-0. 132023E+03 
-0.328320E+03 
-0.328320E+03 

It has been assumed here that all the state variables are 

available at each instant (observability matrix » I,) 

0 

The closed loop dynamics has been simulated and 
Figures (V.3), (V.4) and (V.5) show the transient 
responses to a 6° initial perturbation in roll, pitch and 
yaw, respectively. Figu-e V.3 shows that a 6° 
perturbation in roll is damped out in approximately 13 
seconds. During that single axis maneuver, it should also 
be noticed that the coupling disturbs the yaw degree of 


0.431436E+02 
-0.431436E+02 
0. 132023E+03 
-0 . 132023E+03 
0. 328320E+03 
-0 . 328320E+03 
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freedom, which reaches a maximum amplitude of 0.25° degree. 
Figure V.3a, V.3b, and V.3c show, for the 6° maneuver 
about the roll axis, the forces required from the reflector 
actuators, the efforts produced by the Shuttle's torquers, 
and the components of the equivalent total torque acting 
on the SCOLE system, respectively. The reflector "y" 
actuator and the Shuttle's "x" torquer are the more active 
controllers for this maneuvers, as expected. 

Figures V.4 and V.5 show the response to 6° initial 
perturbation in pitch and yaw, respectively. During the 
maneuver about the pitch axis, the yaw angle is seen to be 
perturbed and reaches a maximum 1° amplitude. The pitch 
angle reaches the same amplitude in disturbance when the 
maneuver about the yaw axis is undertaken. This confirms 
the strong coupling between the pitch and yaw closed-loop 
motions of SCOLE. The pitch maneuver takes about 48 
seconds to stabilize while it takes the yaw maneuver 
almost two minutes to do so. This can be attributed to 
the inertia distribution in this configuration of SCOLE on 
the one hnnd, and to the shorter moment arms (offset 
distances) available to the reflector actuators during a 
maneuver about the yaw axis on the other. Where the 
control efforts are as important, a factor in the selection 
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of a control strategy as the response times are, this 
control law could be chosen over the bang-bang approach.^ 
For the cases shown in Figures V.4 and V.5 the maximum 
control effort from each controller is far below its 
saturation levels. 

▼•2. Control of the Orbiting SCOLE with the First Four 
Modes Included 

This model of the SCOLE is controlled through the 

three torquers on the Shuttle and the six actuators 

located by pairs at z » -L/3; z ■ - 2L on the mast; and at 

3 

Gj » the mass center of the reflector (Figure 1.1). The 
pairs of actuators are arranged in such a manner that one 
acts along the x direction and the other in the y 
direction. The actuators, when activated to provide 
vibration control to the mast, will develop torques about 
the Orbiter center of mass. Each actuator provides a 
maximum of 800 lb^. force; the resulting torque 
contributed by all six actuators is computed as 
-¥■ 

T 1 ' F y L <' , l/ 3+ 2 v 2y /3 + T 3y> 2 ' F x L < v li/ 3 + 
2 t 2* /3 + v 3x ) j - <"« »3x * IF *3,^ 
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This is added to the torques provided by the Shuttle's 

A A A 1 

three torquers: T, ■ M v D, i + M D j + M U k, to 

4 xx 7 7 z z 

yield the total available control torque for the system 
as : 

? - [M x U x + F t L (* ly /3 ♦ 2r 2j /3 + * 3jr )]~i [M y U y - 
F x L(»i x / 3 + 2t 2x /3 + v 3x m ♦ [M, U x ♦ I F y T Jy - T „ 3jt 

After the substitution of the new state variables into the 
generic modal equations. Appendix D, there results: 

•• 

^ - “n + C n - ^ V " ! + - «2 C „ - 3 “o V " 2 


**9 B « ’ “3 V "3 + (a 10 B n ‘ a 4 V " 1 + <«U B n " “5 C n * V "2 


**J2 B n ' *6 C ii + H V "3 + ¥ v 3r s nx ™ * y 2x a <x ^ + T l* V (_L /3>1 

+ V (_L) + v 2 t V «W> * V 1 j V (- W>3 


75 


with the control and global state vectors, respectively, 
chosen as 


U " ^ V lx* v ly * v 2x' v 2y ' v 3x’ v 3y’ U 1 * U 2* U 3^ and 

x " t Y V n 3* Y Y Y V V Y V V Y Y Y t 

vdth » I ItJ !< Is I |v iy | I< 1; and I Juj 1< 1. 

** et t nx » an< i t ny’ n “ ^ * 2, 3, 4, be functions such that 

t nx^ 2 ^ " F x s nx^ z ^ and t nj^ z ^ “ F y s n y( z )» the contro1 

influence matrix can be recast as B ■ 
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t^O-L) 

v» 

0 

0 

0 
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0 

0 

V-W) 

t 4y W3) 

t 4y (-2W) 


V-L) 

0 

0 
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(V. 5) 

The linearized equations of motion describing the 
rotational open-loop dynamics of the orbiting SCOLE 
(modified form of Equations III. 22 - III. 24) are 


"1 - "2 ¥* - v«- “b'V«=-V I = ) - 
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^% n l (I a" V"“o n 3 I «- 3 “5 n 2 r xr + 


I A-^-o 

n-1 


L ‘h V t ^ d 2 n + 


(V. 6 ) 


n 2 I W + A^ H+; ’ 3 % IL - "l “o^ + " 3 “ 0 ^ n - "l 3 ^ 


' 3 « ^JL i U- I i ^,.0 (V. 7 ) 

u n“i n* 1 


n 3 X 2Z " n l r zx " n 2 ^ ^ + “0 n l (I xx ” *JJ + V " n 2 « 0 V 
■* "» “ 0 ^ + 3 "2 -I V - "3 « 0 2 «„ - V + j ^ ** * j 3 K *7 a 


(V. 8 ) 


where , for Che nth node 


d In * M R < Ls ny ( ‘ L) * XY V* L ) + r- f 2 (S n >: 


2n R 0 


M i? u n( Y s nv (-L)); d 


2 , M 


3n “o < r W - W-*)* 
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d 4n - f f l (6 n > + (I K2 + V* 2 + L 2 J)e 0 (-L) 

d 5o " “oV (L9 „<- L > + * •hltX-.yt-L) - 

Ts nx ( * L)): d 7n * M It“o Xr V' L)i and d 8n ■“0 2Ts nx ( * L) 


Equations (V.6), (V.7), and (V.8) can be recast in a 
matrix format as 
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t ? 


T f 


[B,] - tB 1 B 2 J , [Cj] . [ Cl C 2 J 


and 


M m 


0 "W 1, ^ Cr «'‘V I = ) 

■W 0 ajftxr 

•few -w r 0 




*zi *zi ~*x 

d 5I •*& - d 53 -*» 
; d 71 “72 ^ *74 


and, 




m » 

Hw -“ 0 2± e 


‘bi 

"32 

*33 

* 

^4 

- 3 “ 2 r «y “oV 

A * 

N* 

0 

0 

0 

0 

-4i) i 3»_ 2 1 - uj 2 n: -i 

OB 0 yz 0 u xx yyj 


> 
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From the generic modal equations, the expressions for 
•• 

the different A^, (open-loop) can be expressed as follows: 


A 1 


'* 1 ' 




CN 

5 «< 


A 2 
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r ..1 


■H 


+ N 


+ 

•i 

C 2 

A 3 


A 3 


n 2 


A 4 


v 

**• 

, n 3 . 
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with 



where u» i is the frequency of vibration of the ith mode, 
and 
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After substitution for the open-loop A (i»l - 4), 
Equation (V.9) can be rewritten in the following matrix 
"format : 


[ A ll 




"1 

n i 

tt 

■ - 
• 

n i 


s 

+ f A 2 C 2l 

n 2 

n 

+ t A 2 C 2l 

• 

n 2 

• 

+ t B l] 

”2 

• 

n 

1 

3 J 

* 


n 3 

» * 


. 3j 


; - 





h 

• 

V 

V 

V 

A l‘ 

CM 

< 

+ [Cj] 

n 

2 

+ £C 2 ] 

A 2 

•*3 


n 

3 


A 3 

A 4. 


» — 


A 4 


(V.10) 


which, when the global state vector is taken as: X^* » 

[Hi* n 2 * Hgt A^ , A2 » A^ # A^, n ^ , r» 2 * ^ 2 » A^» A^], 


is equivalent to 


[- Cj ) I -k\(^ 2^)1 - b' 2 ] X 


.-1 


- 1 / 


-u* 


The new state matrix A, for the SCOLE system 


therefore, is 
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[ 0 ] 

[ 0 ] 


[ 0 ] 


[ 0 ] 


-A 1 "* 1 (A 2 C 2 +C 1 * ) | -A 1 “ 1 (A 2 B 2 +C 2 ') 


3x3 


[c 2 ] 


3x4 


4x3 

[S 2 ] 

4x4 


[13] j[0] 

I 

J_ 

[0] : [I4] 

—t * • * i i » 

-*1 (*2 C 2 +B l > - A 1 B 2 


3x3 


» * 
[C ] 
2 


4x3 


3x4 I 4x4 

Here again, the control U which minimizes the performance 
index 




o* 

5 


T T 

(X QX + U RU) dt is obtained after using the 


ORACLS package to solve the steady state Riccati matrix 

9 

equation* Figure V.2 shows the same type of parametric 
study previously conducted for the model of the rigidized 
SC0LE. Since it is anticipated here that large amplitude 
slew maneuvers of the flexible SC0LE will be conducted 
based on the control strategy developed herein , the 
concept of split weighting of both the state and control 
penalty elements was considered. The criteria of 
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selection being based on the control strategy capability 
to slew the flexible SCOLE through large amplitude angles, 
^ ® » and $ equal to 20° respectively without any 

of the actuator reaching saturation level. 

It should be noted here that the state and control 
influence matrices have dimensionality of (14 x 14) and 
(14 x 9), respectively. 

The equations describing the closed-loop system, 

X ■ AX +BU have been numerically integrated and the 
corresponding mathematical model simulated for 

<Hiag.[5xlG 6 ,5xl0 6 r 5tl0 6 ,5xl0 4 ,5xl0 4 ,5xl0 4 ,5xl0 4 , 10,10, 10,10,10, 10, 10] 


and R as diag [10,10,10,10,10,10,1,1,1] 

The transient responses to some initial 
perturbations, depicted in Figures (V.6) - (V.17) confirm 
the controllability of the flexible SCOLE system. During 
the simulation of this model, the three attitude angles 
(roll, pitch, and yaw) are each subjected to a 6° single 
axis maneuver. For each case, the effects of such 
displacements on the modal amplitudes of the first four 
modes are shown (Figures (V.6) through (V.ll). 
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The largest disturbance in the flexible modes, caused 
by the variation of an attitude angle is observed during 
the roll axis maneuver, the first mode is the most 
excited; its amplitude doesn't exceed 0.13 ft. (0.1% of 
L) . All the transients are damped out within 25 seconds. 

This is due to the contribution of the additional 2 
pairs of actuators located on the mast at - -L/3 and z 
* — 2L/3. During the 6° maneuver (from equilibrium) about 
the roll axis , the reflector "y" axis actuator provides a 
maximum of 210 lb. while the forces in the two "y" 
actuators located a z ■ L/3 and z • -L/3, ‘reach 120 to 

80 lb, respectively. The Shuttle "x" torquer provides a 
maximum of 2800 ft. lb torque bringing to 52,500 ft. lb 
the maximum value of the x component of the composite 
control torque required for this maneuver. This, when 
compared with the total maximum torque of (35,000 ft. lb) 
required during the same maneuver of the rigidized model 
of SC0LE (figure V.3c), shows an increase in the total 
control torque of 50%. However the reflector "y" 
actuator, when flexibility is included, provides less of a 
contribution than for the rigidized case. 

It should be noted here that because of the 
additional pairs of actuators located on the mast, one can 
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now exploit the moment arras provided by the two reflector 
actuators without the same risks of perturbing the pitch 
or roll which exists for the rigidized model. The 
actuators located on the mast would help to prevent such 
undesired displacements in pitch and roll. 

In turn, each of the four flexible modes were given 
an initial amplitude equal 1.0% of L, to stay within the 
linear range. Figures (V. 13) , (V. 15) , (V.17) and (V.19) 
show the transient responses to those initial 
displacements. Also depicted is the result of intra- 
flexible modal coupling. For this control strategy, the 
disturbances in each of the flexible modes, for the 
initial conditions considered herein, are damped in 15 
seconds while their effects on the attitude angles take 
almost 25 seconds to disappear (see Figure (V.12, (V.14), 
(V.16), and (V.18)). The effects of the coupling between 
the flexible and rigid rotational modes is best observed 
in Figure (V. 16) (roll response) when the system is 
initially excited in- its third mode. The control effort 
required here is similar to that of Figures (V.6—V.11) 
described above. 
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Preliminar2j__Slev_Jlajieuver^ 

In this section, the equations governing the motion 
of the rigidized SCOLE, outside of the linear range, are 
developed from the most general rotational equations of 
motion previously derived. The control laws obtained from 
the application of the linear regulator theory to the 
linearized model of the rigidized SCOLE are tested for 
^®rge amplitude manuevers. The closed loop system 
dynamics are numerically simulated. For single axis slew 
maneuvers about the roll, pitch, and yaw axes, 
respectively, the time responses- for the Euler angles, the 
control efforts required of the reflector actuators, 
control torques demanded from the Shuttle's 
torquers, and the components of the total control moments, 
are depicted in the subsequent figures. This enables one 
to determine the margin left in which to optimize the 
control strategy without causing saturation of the 
controllers. 

In the absence of flexibility in the system, Equation 
(III. 18 ) becomes: 

H syst/G * & x^l ~ Q z l SA + Q xj M + V** + cY) + ^1%^ 
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+ { fl T ^2 + If. 2 a 7 - f^(at + dD + a ^ j 

+ (8 z % 3 -n r % 4 + H t CK-ar)+ n^k (V.U) 

where ft^, ft y , and ft^ have been defined in Equation 
(III. 3). In the absence of flexibility, [T R R2 ] - [I 3 ] 
which is equivalent to - T 22 - T 33 - 1.0 and T\ ^ m 0, 

it j • 

Under these new assumptions 

fl l * R x ; fi 2 " n y ; and Q 3 * z 

a «• ii L - «y 

7 2 

b ■ Q „L + U and 
X z 

c ■ ft Y - ft X 

x 7 

Therefore, Equation (V.ll) can be rewritten as 

Vo ’ { a x¥Vf 2 + W? * *»■ a r ^ H 

+ {a y a S2 +£0, 2 + i R2 + m,j(l 2 +x 2 ))- a x ^ rr + a z 

+ (ft z (I S3 +I R3 +M r (X 2 +Y 2 ) + ft x M R XL + ft y M R YL}k - 0 


(V.12) 



88 


which can again be recast as: 

H »fse^G * "x 1 + h t J + H z k (V.13) 

The equations governing the motion of the rigidized 
SCOLE system during large amplitude maneuvers in the 
presence of gravity gradient and control torques are* 
obtained as: 


i) The Roll Equation 


• 


H * + Ay H z " «z H y " T x 

(V. 14) 

ii) The Pitch Equation 

K * a, H x * a, "z - 

(V. 15) 

iii) The Yaw Equation 


■« + H y * H x ' T z 

(V. 16) 


where T^, T^, and T z are the components of the external 
torques acting on the system (including the control 
torques previously derived for the linear model of the 
rigidized SCOLE where the feedback now depends on the 
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Euler angles and their rates for maneuvers made 
relative to the Shuttle roll, pitch, and yaw axes). 

The closed-loop system dynamics described by 
Equations (V.14), (V.15), and (V.16) have been numerically 
simulated and the results are shown in Figures (V.20) t.o 
(V.32) . 

Figures (V.20), (V. 25), and (V.29) show the time 
responses to an initial 20° alignment in roll, pitch, and 
yaw degrees of freedom, respectively. It is seen that a 
20 slew about the roll axis can be achieved in about 30 
seconds. The same amplitude maneuver about the pitch and 
yaw take 45 and 100 seconds, respectively. This is due 
to: 

1. the system inertia distribution 

2. with equal amounts of torques available in the 
Shuttle for each maneuver, the roll and yaw 
maneuvers benefit more from the actuators 
located on the reflector for which the length of 
the beam is then a moment arm. 

For this control strategy, each of the single axis 
slew maneuvers about the roll and pitch axes used 85-90% 
of the control forces available from the corresponding 
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actuator located on the reflector, and 80% and 60%, 
respectively, of the control torque available from the 
corresponding Shuttle torquer. None of the controllers 
reach saturation. The control strategy which has been 
designed so as to avoid an excessive use of the actuator 
forces, relies mainly on the Shuttle's "z" torquer to slew 
about the yaw axis. It is seen for this maneuver (Figure 
V.31), that the corresponding Shuttle torquer is used at 
99% of its maximum capacity. Also depicted in the Figures 
(V.24), (V.28) and (V.32) are the components of the total 
control moments for each case (moments of the reflector 
control forces taken about the Shuttle's mass center, plus 
moments of the Shuttle's torquers). This will make 
possible a comparison between this strategy and other 
future control laws which would be based on the two point 
boundary-value problem, where this or combinations of 
control inputs may be employed. 

In conclusion, it is seen that a control strategy 
derived from the linearized model of the rigidized SCOLE, 
based on the linear regulator theory, works well when it 
is used for single axis slew maneuvers through amplitude 
angles as large as 20°. 
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figure V. 17 LINEAR MODEL OF SCOLE WITH FLEXIBILITY 
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CHAPTER SIX 

CONCLUSIONS AND RECOMMENDATIONS FOR FUTURE INVESTIGATIONS 

In this thesis, it has been seen that: 

1. The SCOLE system, with gravity-gradient torques 
included in its open— loop dynamics, is unstable* 
This is due to the inertia distribution of the 
system in the particular configuration where the 
Shuttle roll axis nominally follows the orbit. 

2. The equations describing the pitch motion 
decouples, within the linear range, from the roll 
and yaw equations, when the gravity-gradient 
torques effects are present in the system 
dynamics, and when the system is without offset 
or when the offset is parallel to the roll axis. 

A result similar to the one depicted here was 
derived for the tethered platform system. (11) 

3. In the absence of control forces and torques, the 
system will oscillate about a new equilibrium 
position , The amplitudes of the oscillations 
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grow with the offset distance and the frequency 
of oscillation during that motion depends on the 
frequencies of the modes taken into account in 
the model. 

4. A control law judiciously derived for the linear 
model of the rigidized SCOLE can be used to slew 
the system without reaching the saturation level 
of the controllers. It is anticipated that the 
trade-off between maneuver time and overall 
control effort would be in favor of such a 
control law as compared with the bang-bartg 
strategy or the two point boundary value problem 
approach. 

The coupling between the elastic displacements 
and the rigid modes is strong enough to suggest 
more accuracy in modeling this class of offseted 
and large flexible structures. 

The author suggests the following topics for future 
research. 

1. In the case of the rigidized SCOLE model for 
single offset parallel to the .-oil axis, the 
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equation describing the pitch motion of the rigid 
SCOLE, decouples from the roll and yaw motions. 

In such case, a control law could be derived 
analytically and compared with the control law 
derived using the linear regulator theory when a 
maneuver is done about the pitch axis. 

2. Since the SCOLE design* challenge consists of 

slewing the SCOLE configuration in a minimum time 
through a 20° line of sight angle, it is 
conceivable that the high rates at which the 
slewing maneuvers occur would modify the 
structural configuration of the system, at least 
for those short periods of time the maneuvers 
would take. Therefore, i study could be 
conducted on a model which would include a load 
equivalent to the effect of such induced 
centrifugal forces and the subsequent modal 
shapes and frequencies could be compared with 
those at hand at the present time. If signifi- 
cant differences in mode, snapes/f requenc" 1 es 
during slew maneuvers are observed, then the 
slewing simulations reported here for the 
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rigidized model should be repeated and compared 
both with the uncorrected and also with the 
corrected flexib le models, 

3. A study could be undertaken which would derive a 
global control law compatible with the two-stage 
strategy during which first the system would be 
slewed as if rigid and second the induced elastic 
vibrations suppressed. 
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STRUCTURAL ANALYSIS OF THE SCOLE SYSTEM 


In this Appendix, the mode shape, the corresponding 
frequencies, and modal amplitude are derived for the SCOLE 
system. In this analysis, it is assumed that the beam 
mast has: 

a) a uniform density; 

b) a circular cross section; 

c) a uniform distribution of stiffness; 

also that the displacements and slopes are small. 


GoTarning Differential Equations 

The governing partial differential equations for the 
beam are comprised of two one plane bending equations. 


(A.l and A. 2), and one axial torsion -equation, (A. 3). 
For the x— z plane bending, one has 


3 1 


3 a ( z . t) m El 
2 p A 


u (is , t) 


(A.l) 


where p is the density of the beam, A its cross sectional 
area and El its x-z plane bending stiffness which is 
assumed , in what follows, to be equal to its y-z plane 
bending stiffness (circular shaft) 
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Assuming the following form for u(z,t.) (separation of 
variables) : 

- P x (t)s x (z), 


equation (A.l) can be rewritten as 


IK («) 

(4) 


or 


El _jc_ 
pA 3 


This equation is true if, and only if both sides are 
equal to a constant, say, - u 2 , yielding 


P x + <u x P x - 0 


which integrates into: 


P x (t)« cos(uj x t+a) where a is a phase angle, 
right side, 

( L ^ 7 n A 

0 


From the 


s < 4 > - w 2 £A s 

x x El x 


L.tclng 

form: 


, the general solution has the 


3 x * ^1 sin6 x z + cqs8 x z + C x sinh8 x z + cosh8 z 


u(z,t) - cos ('o x t + at) {A 1 sLi 3 x z + Bj cos8 x z + CjSinh B^z + ^ coshB^z} 
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The y-z plane bending Is described by: 
i 2 FT a 4 

2 (v(z.t)) - ff -2-r- (v(z , t) ) 

3t pA 3 Z 


Assuming v(z,t) of. the form 
- P y (t) s y (z) 

After substitution , Equation (A. 2) becomes 


(A. 2) 


-s 


7 




P 


y 


» ei 
" pX p y 

El a v ( 4) 
P A g 

y 


s (4) 

y 


or 


which is true only if both 

sides are equal to a constant say, - w y* A reasoning 

similar to the one used to analyze the x-z plane bending, 
yields 


v(z,t) - oos(w y t+Y) {AjSlnS 'z + cos S^z + sdiihBy z + D 2 co *8 

where 8 4 » a> 2 

y Ely 

Finally, the z axis torsional bending is described by 


(A. 3) 
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where 6 is the modulus of rigidity of the beam. Here 

again* assuming the separation of the variables possible: 

♦ (z,t) P z (t) e(z), and letting g' 2 « p / G , there 

z 3 

results: 


* (*,t) 


cos (w 2 1 + 5 ) 



f 

sin 8 z z 


t 

+ B. cos 6 z} 
J z 


The equations giving u (z,t), v(z,t), and $ (z,t) are 

•more convenient to use when the position variable is 
transformed into a nondimensional form. For this reason, 
the variable, e" z/L, where L is the length of the 
undeformed beam, is used.- After substitution into 

v(z, t) , and *(z,t), those equations become: 


u(e,t)-cos (o) x t+a) {AjSinB x e + BjCosS x e+ CjSinh 8 x £+ D^osh 8 x e} 
v(e ,t)-cos( aj yt+y) {AjSin 8 y e+ J^cosgy e+ C 2 sinh 6^6+ D 2 COSI 1 B y e } 

* (£,t)«cos (<* t + <5 ) { A~sin 3 e + B-cosS e} 

Z Z j z 


m -P A ■ ■ Ul 2 
x (EI) x X L 


PA u> 2 L 4 


y (EI) y y 

« 2 - p - 2 T 2 

p_ “ 7T to L 
Z G 2 


with. 
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In this investigation it is assumed that 


<*> z and P x (t) - p y (t) - p z (t) - p(t) 


6 - - 6y 
8 - if El 

z L GA 


GA 

EZ 


S z L 


or 


0) a 


X 


0) a 

7 


again 


Boqndarr Conditions 

In our model, the offset of the mast attachment point 
from the center of mass of the reflector, along with its 


products of inertia. causes a kinematic coupling between 
the displacements in the different degrees of freedom. 

The following relationships between shear, moment, and 
beam displacement are used in the boundary conditions ^ 

V a- ~ (e,t) - - El i 2 £< 2 ’ t ) 

x 3e 3 3z 3 


V 

7 


_ El ( e , t) 


L 3 3 e 3 



(z , t) 
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where: 

V x » shear force in the x direction 

V m shear force in the y direction 
7 

M^, My, and M z • moment components about the x, y, 

and z axes, respectively. 

I ■ polar moment of inertia of the beam. 

P 

Let Mg be the mass of the Shuttle and M R that of the 
reflector. The shear force at an end of the beam is 
assumed equal to the mass of the corresponding body at 
that end multiplied by the acceleration of that end; and 
if we also consider that the displacement in the x 
direction of a point located at z*0 is given by 
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u(0,t)-DyQ $ (0, t) and that in the y direction by 
▼ (0,t)+DxQ (0, t) , (Dx 0 and Dy Q are the displacement 

components of the centroid of the cross sectional area of 
the beam at z»0), therefore, at the Shuttle end, 

V x l- M s p {s x (0) “ Dy o 9 (0)} “ " ^ s x 3) <0)P(t) 


V 7*"n M a p{S y (0) " Dx 0 (0) m - ^ 5 C3) (0) P(t) 

7 e»0 7 l J y 

Taking into consideration the fact that p - p(t) 
and substituting it into the shear equations yields 

- u 2 M s [ a x (0)-Dy 0 9(0)] - s ( 3> (0) 

L 

- 2 M s t» y (0)-Dx 0 (0)] - -Si 

L 


but 


0) 


6 


4 El 


pAL 


After rearranging terms, one obtains 
•xVorrfL 8 4 [s x (0)-Dy 0 l,;0)i 
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(3) M « 4 

S y (0) " pAL 6 Dx 0 9(0)] 


A similar reasoning at the reflector end of the beam 
taking into account the equilibrium of the beam would give 
(3 > M 

®x (1) - ^AL 8 t- » X U) + Dy 9(1)] 


.(3) 


M. 


*T (1) ’TaL 8 t- > Cl)-DX L 9(1)1 


Bending Momenta 

The next four boundary conditions involve the 
moments, M xq , M yQ , , and M yL _ at z - 0 and z - L, 
respectively, on the beam. Assuming the nonlinear 
coupling and all the products of inertia, except I , to 

X z 

be negligible, one may write 



• # 


xx 


+ 





M » I e + i e 
y yy y *y x 


where and 0 y are the angular accelerations of a point 
on the mast about the x and y axes, respectively. 
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since 




M - El u (2) 

7 t 2 


from Equations (A) f 


and the angular displacements e , and e are given as 

X 7 


0 

x 


_1 . 3 v 
L 3e ; 


and & 

y 


1 3y_ 
L 3 e 


substitution of the general expressions for u(e ,t) ,v(e, t) , 
v($,t) and (e,t) into the boundary conditions (A and B) 
at e ■ 0 and e * 1, respectively yields 


s 





s< 1 > 

y 


( 0 )) 


si 2) C0) 


pAL 


J(-I S2 


(It should be recalled that the term I as applied to the 

xy 

Shuttle is zero.) 


( 2 ) 

s y (1) 




(1) 

S y ~ M^DxDy 



( 1 )] 
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( 2 ) 


s_ 


( 1 ) 


pAL 


3 ^ X R2 




DxD 


S (1) l 

y 


vhdr6 I s ^ | 1^ ^ i Xg 2 1 stnd 2 th6 moi&6iits of 

inertia, about the x and y axes, of the Shuttle and the 
reflector, respectively. The x and y axes considered here 
pass through the respective interface points. 


Torsional Moments 

These moments are caused by the masses and moments of 
inertia, about the z axis, of the end bodies: $ (z,t) I 

zz 

They are countered by the beam internal moment given by 


M - GI 3$ 

2 p 77 


Writing the equality between these two moments yields: 

at the Shuttle end 

B 2 

9(1)(0 > • St!" 'S3 9 (0> <?> GI p If ■ 


~t the reflectot end 

.2 

. * 
e 


ri'i b" 

^ •’(1) - _5_ [-X r39 (1)+M r Dx l 3 <1)-M r D 7l 3 x (I)] 


pLI 
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I 


where Ig 3 and Ijj 3 are the moments of Inertia of the 

Shuttle and the reflector, respectively, about the z axis. 

Substituting the expression for s (e ) , s (e ) , and 0( e), 

x y 

and the appropriate values of e at the boundaries into 

the boundary equations, the following ten linear 
equations are obtained. It should also be noted here that 
" DJq ■ 0 since at the Shuttle end the centroid of the 
beam cross section coincides with the center of mass of 


the Shuttle, 


-A 8M s 

A 1 * ( — 2 .) 
1 l pAL ; 


SM 


B 1 + Cj -(_£) D, . 0 


PAL 


(A. 4) 


B 2 + C 2 D 2 - 0 


(A. 5) 


A,( 


% 

pAL 


sing- cos 6} + B^{ 


PAL 


cos + sin } + C. 


gM_ 

{L* 

• pAL 


sinhS-coshB } 


BM 


[ R eM n 0M_ 

X cosh0fsinh0} + A 3 { sin0 2 } + B^ Dy^ cosB^ 

(A. 6) 
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• 8M gM„ 

A»{ — r-^sin - oos g} + a,{ — - 

* pAL p pAL 


cos.ef 




PAL. 


8M, 


co*B«tah6) + *3(i^S “"L « to l> + %( °«L «• V - 0 W 


*•> C- 8 — + C 2 (- |“ ) +D, -. 0 


PAL' 


PAL' 


CA.8) 


B 3 ! 


Al < 17 ? > - b i + 0 i ( ^ +D i- ( ^r 'W ■ « 


8 3 I 


S2, 


81 


SI 


(A. 9) 


(8 3 M r Dx_Dx) 
pAL 3 


(AjCOsS-BjSinS+CjCoshB+DjSinha) 



sin8-cos8 } 


-C 


2 


(— ,Ig,eo»he+sitth6) - D,(J2 i d . 
8 AL 3 3 pal 3 R1 


sinhg+coshe } 


0 (A. 10) 


-Aj{JL.Ir 2 cosB+sinS} + sin8-cos8) 

PAL J PAL 3 

+Cj {sinh6-_8fl R2 coshS) + Dj {coahB- ^ 3 I R2 sinhS} 
PAL 3 pAL 3 
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+ (g^M ^Dx Dy ) (A 2 cos S-B 2 sin ^+C 2 coshe+D 2 sinh 6) » 0 (A. 11) 
pAL 3 


A 3 + B 3 " 0 

P 


(A. 12) 


(Sj Dy M R ) { Arsing +B^cos8+C^sinh8+ D^coshS } 

pLIp 


»(8 t Dx M r ) {-A 2 sin6+ B 2 co36- C 2 sinh6 -D^oshe } 

pLI 

P 

+ A 3 1 6^3 8in6 z + cos8 z ) - B 3 {8 z I R3 cos6 2 +sin6 } - 0 (A. 13) 

PLI P 

Equations (A. 4) through (A. 13) can be recast in the 
following matrix format 

[Z(B)]{[A 1# B^ C 1# D r A 2 , B 2 , C 2 , D 2 , A 3 , B 3 ] T ) -[0] 

where [Z(6)J is a 10 x 10 matrix whose entries are 
functions of 6. Mon zero solutions for A^ B i , C^, and D ± 
exist only when the determinant of [Z(8)] is zero. 

The equations derived here being identical to those 
derived in reference 3, the values of A ± , B i , C^, and D i 
obtained therein have been used (table A.l) and the 
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projections of the first four mode shapes plotted: Figure 

(A.l) through Figure (A. 12). 



•hM* k 1 


6 























































MODE SHAPE X-Z PLANE BENOING 
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e = z/L 

FIGURE A-ls PROJECTION OF THE FIRST MODE SHAPE ONTO THE X-Z PLANE 
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3 

i 


e - z/L 

Figure A-2: PROJECTION OF THE FIRST MODE SHAPE ONTO THE Y-Z Plane 





DIANS) 


MODE SHAPE X-Z PLANE BENDING 
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FIGURE A- 4 j PROJECTION OF THE SECOND MODE SHAPE ONTO THE X-Z PLANE 



MODE SHAPE Y-Z PLANE BENDING 
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FIGURE A- 5 ; PROJECTION OF THE SECOND MODE SHAPE ONTO THE Y-Z PLANE 
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MODE SHAPE Y-Z PLANE BENDING 
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FIGURE A- 8 : PROJECTION OF THE THIRD MODE SHAPE ONTO THE Y-Z PLANE 
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i 


FIGURE A-9: TORSIONAL DEFLECTION OF THE THIRD MODE SHAPE (0 IN RADIANS) 


MODE SHAPE X-Z PLANE BENDING 
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FIGURE A- 10 : PROJECTION OF *iHE FOURTH MODE SHAPE ONTO THE X-Z PLANE 
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ORIGINAL PAGE IS 
OF POOR QUALITY 



FIGURE A-Il : PROJECTION OF THE FOURTH MODE SHAPE ONTO THE Y-Z PLANE 




TORSIONAL MODE SHAPE 
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FIGURE A-12: TORSIONAL DEFLECTION OF THE FOURTH MODE SHAPE (0 IN RADIANS) 
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APPENDIX B 

ANGULAR MOMENTUM OF A RIGID BODY 
Transfer Theorem ^ 

Let B be a rigid body with its center of mass located 
at point G and let I be some arbitrary point. 



By definition, the angular momentum of B about the 
point I is given by: 


t 


I 



x $(P)dm 


(B . 1 ) 


where V (P) is the inertial velocity of the dm located at 
point P. 


The vector IP can be rewritten as 
fp - fG + GP 


(B. 2) 


and 


V(P)«V(G) + ax &> 


(B.3) 



I 


160 


where V (G) is the inertial velocity of G an<L ft is the 
inertial angular velocity of B. 


IP x V(P) can therefore be expanded as 

x f( P) - & x \f(G) + jfex(x$)+&x V(G) + GP x <&x GP) (B.4) 

Each of the terms can be integrated as follows: 




IG x V(G)dm - IG x V(G) J dm - MIG x V(G) 

M 


f 


M 

^GP x V(G)dm - -V(G) x JGPdm - 0 

M M 

for, G i.a the mass center of the body B. 


■ 5 - 


Jig X (^xGP)dm • (IG x?) xf GP 

M 

because of the reason stated above 


Let now GP '« Xi + Yj + Zk 


V + V + V 


a * 


(B.5) 


(B.6) 


(B . 7) 


(B.8) 


(B.9) 


where i, j, k are any convenient reference axes fixed to 
B. 


GP x (^xGP) - | GP | - (ff.GP)GP 


(B.10) 
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Equation (B.10) after substitution of Equations (B.8) 
and (B.9) becomes: 

GP x («xGP) - [(Y 2 +Z 2 )A - XY® - XZ 0 ] i 

x y z 

+ [-XY« + (X 2 +Z 2 )fl -ZYA ]j 

x y z J 

+ c-zxn x - YZfi y + (X 2 + Y 2 )n z ]k (B.ll) 


Because 


f 


(X 2 +Y 2 )dm - I 


ZZ/G 


M 


1 


XZdm ■ — X 


XY/G 


and 


r* 

XZdm - -I ZX/G 
M 


J 


M 


YZdm - -I 


YZ/G 


I 

f 


(X 2 + Z 2 )dm - I 


YY/G 


(Y 2 + Z 2 )dm » I 


XX/G 
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it can now be seen that 


GP x (fix 3P) dn ■ I G n (B. 12) 

which is the angular momentum of B about G. Gathering the 
terms yields 




H G + M IG x V(G) 


(B. 13) 
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APPENDIX C 

GRAVITY GRADIENT TORQUE ON THE SCOLE SYSTEM 

In what follows here, the gravity gradient torques 
will be derived for three different configurations of the 
SCOLE system. In the development of the expression of the 
gravity gradient torques for all three cases, it will be 
assumed, without great loss of accuracy, that the center 
of mass of the entire system coincides with that of the 
Shuttle Orbiter. Under that assumption, the unit vector, 
a, defining the local vertical can be expressed in the 
Shuttle body frame as the following functions of Euler’s 
angles (See Chapter III, Section 1). 

A A 

a » sinQcos 9i - (cos Qsin <|;+sine sin<j> cosiji ' j + (sinQ sin<)> sim|> - 

a 

cosecosij/ )k (C.l) 

Expression for the Gravity Gradient Torque. N 

* - 3 “0 • x I syst/G' ® < C - 2 > 

where, u>q is the Shuttle (circular) orbital angular 

m 

velocity; I S y at /Q# the inertia tensor of the system at its 
center of mass. 
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The total tensor of inertia about G is the sum of the 
inertia tensors of the components transferred at G using 
the parallel axis theorem. 
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SCOLB System Without Offset 


'syst/G 


-I, 


where 


U - 


0 

I. 


hi + X R1 + -Mk 2 + M r L 2 

l S2 + X R2 + -yi- 2 + M R l2 

r S3 + *R3 and *4 * *S4 


Under Che esaumpclon of small angle approximation for 

the Enler angles, there results 

-*• 

N - 3a» 2 0 ( * (I 3 -I 2 )i+[-I 4 -e(I 1 -I 3 )]j+I 4 '<-k} (C.3) 

SCOLB System with Offset in the "I" Direction 


‘syst/G 


0 

I 


-I, 


-Ii 


with 


V 1 


4* W L : I 2 - X 2 +M R X " 
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and 


-i 3 +m r x 2 


again 


» -3«? 0 t*(I 3 ' -I 2 ’ 
(C.4) 


)i+(-I 4 ' - e(I 1 ’ -I 3 '))j+I 4 ' * k) 


Actual SCOLE Configuration in the Pndeformed State 


In the actual configuration, offset in both the "X" 
w y w directions, the design challenge paper^^ provided 
the Inertia tensor of the whole system as 


xx 


-I 


yx 


-1. 


xz 


‘syst/G 


-I 


xy 


-I. 


xz 


which yields 


yy 

-1 


yz 


-1 


yz 

I 


zz 


N - 3u {-I ^ + d»(I -I )+ ei ]i 
jz zz yy' xy J 

♦ t-x„+*i XT -9(x„-r„)]j+(-9i_ + M„)ic) (c.5) 


XX zz 


yz xz - 
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APPENDIX D 


Generic Mode Equations 

Consider an elemental mass, dm, of the body whose 
Instantaneous position from the center of mass of the 

* 4 * 

Shuttle is r. The equations of motion of dm can be 
written as 


adm - L(q) + fdm + 7dm (D.l) 

where a is the inertial acceleration of dm; 7, the 
gravitational force per unit mass; 7, the external force 
per unit mass; q, the elastic displacement of dm; and L, 
a linear operator which, when applied to the small elastic 
displacement q, yields the elastic forces acting on dm. 

The gravitational force per unit mass, 7, can be 
expressed as 

7 - 7 0 + M? (D. 2) 

where f Q is the gravitational force per unit mass at the 
center of mass of the body considered and M, a matrix 
operator. 
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In what follows, the generic mode equations will be 
derived based on a Newton-Euler formulation. The 
principal assumptions made in this development are; 

1. within each component of the system, the mass and 
structural properties are uniformly distributed; 

2. the material of each component is isotropic; 

3. the system is deformed in such a manner that it 
experiences only small strains (within the linear range). 

4. the elastic displacements are small as compared 
with the characteristic linear dimensions of the system; 

5. the natural mode shapes of free vibrations of the 
system are known a priori; 

6. the system is nominally earth pointing; 

7. the system is considered to be closed; no mass 
transfer across its boundaries. 

The vector equation (D.l) can be rewritten in the 
frame moving with each body as: 

— * * 

£ a ou + r + + + ux(wxr)] dm » L (q) + (7+e) dm (D.3) 

• • 

Note that r and r are the velocity and acceleration of dm 


i 
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as seen from the body fixed frame. The symbol cu refers to 
the inertial angular velocity .of the body. The instantan- 
eous position vector, r, of dm can be written as 

^ » r o + ? (D.4) 

where r o is the position of vector of dm with respect to 
G, center of mass of the Shuttle, in the undeforraed state; 

and q is the elastic displacement of dm. Hence, 

• • •• •• 

• r • q and 7 - q* (D.5) 

For small amplitude elastic displacements, one can write q 
as a superposition of the various modal contributions 
according to 

OB 

T - I A n (t) \ <7 0 ) CD. 6) 

n»i 

where A n (t) " P n (t) » modal amplitude 

- an* ♦„(«•„) - S XB i * S 7n j + 8 n k CD. 7) 

The mode shape 4> n (r Q ) is associated with the 
natural frequency, u> n , and satisfies the following 
conditions : 
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♦ 


m ’ 


4 dm 
n 


6 M 
on n 


(D.8) 


where M n is the generalized mass in the nth mode. 


and 


2 — 

L( <b ) « - u <b dm 
n' n n 


i 

s 

M 


<P _ dm * 0 
n 


r x y dm 
o n 


(D.9) 
(D. 10) 

(D. 11) 


This here assumes that the fundamental structural 
frequency, , is much greater than the orbital angular 
velocity, ■ 0.0011 rad/s, and enables one to use, with 
a high degree of accuracy, the mode shape functions 
corresponding to a non-rotating structure. 


Generic Mode Equations 

The generic mode equations are obtained by taking the 
modal components of all internal, external and inertial 
forces acting on the system, i.e.. 


$ 


n 



*• • • 

[a* + r + 2i xF + "3>xr 

cm 


• [LCq")/^ + f + e] dm 


+ “> x(^xr)l im ■ 


(D. 12) 
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The various terras appearing in equation (D.12) can 
now be expanded as follows: 


<b . a dm * a . 

n cm cm 




dm - o (D. 13) 


1 
M 

|? n . r<ta - f ? n . f <ta . J. *„Ct)T m ) dm 


See Eq. (D.10). M 


M M . M 

because of the result established in Equation (D.8), 

• • •• 

. r dm . A n M„ CD- 14) 

. (2TJxr) dm - 2 fT n . (uxf) dm (D.15) 


. M 


j*7 n . ( wx 


n . (wxr) dm - j *n * ( wxr 0 ) <*■ + • faxq) dm (D.16) 

M M 


f ♦ n • «x(5uxr) dm ■ ^ d» n . u)x(u>xr o ) dm + . wxCwxq) 

M M Af 


K • 


1* (q) / j_ dm ■ — ui A M 
/dm n n n 


dm (D.17) 


(D. 18) 


M 


I 


T • fdm 

a 


r n • 

^ uk 


edm ■ E 


^♦ n • V® + n * M< F 0 )dm + ^n * W*" (D ’ 19) 
M M M 


(D.20) 


where E^ is the modal contribution of the external 
forces (control forces) in the nth mode. 
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After substitution of the values for the. integrals 
into Equation (D.12) and rearrangement of the terms, the 
generic mode equations are obtained in che following form 


A + % A n + * /M + l 

n n n n n i.i 


* /M 
mn n 


t8. 


I 

m« 1 


8 


mn 


+ .E 


D J /M n 


(D.21) 


where 




xr Q ) + <ux(«nxr o )] dm 


M 


y * r- r -- -i- - — 

" ^n * t^wxq + tox< l + “•xC^xq] dm 


I * . fi 

«-i “ J 

M 

g n * (*n - M 

L *»" *J \ ,M <’ 


(r ) dm 
o 


)dm 


(D. 22) 

(D.23) 

(D. 24) 
(D. 25) 



(D. 26) 
(D. 27) 


M 

Because, in Appendix A, the SCOLE system was assumed 
to be an unconstrained structure with end masses having 
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inertia, the orthogonality conditions expressed in 
Equations (D.10) and (D.ll) result in: 


D n * -j"n.y.t ? i.- f 0 d “ ' 

- r 0 - Jflayat * n dn * ° 

It is assumed here that the transverse displacements 
are small as compared with the characteristic dimensions 
of the system. In this first approximation analysis, 
t®rms involving the integral of <fr.q will be assumed small 
as compared with terms involving $ .Tq and, thus, 
neglected. As a consequence of this 


ao 

i 

m»l 


$ 

ran 


leaving 


00 

l 

m-1 


’ran 


0 


(D.28, 


A +o) A + 
n n n 


V n 


(8„ + E n )/M n (D. 29) 


Since the control forces consist of six actuators 
locr_ted in pairs at 6 ^ , the reflector ma^s center, and at 
two points on the mast corresponding to z » -43.3 ft 'and z 


-86.6 ft, 
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E n " M » F x (s xn ( - L > + 3 xn <’ L ,3) * 3 xn (- 2L '' 3 )> 
* Vr (S Tn ( - L) + 3 7 „<- L /3 > + 3 yn <-2L/3>) 



0 


“ tl3r ' 3 xn * A ln aia V + B l„ coaB „ 3 + c ln 3 *»*> »„* + D ln CortB n x 

3 r» " A 2n 3ln6 n 2 + B 2n c036 n z + C 2n 3lnh V + D 2n c03h8 I1 2 

» » 

and ® n ■ Aj n aln 6 z + cos 6 z 

are the x, j- components of the mode shape vector and the angular 
displacement of a point on the mast about the z axis, respectively 

’* lth » - 8 nfl? 

♦ can be rewritten as 
n 
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*n " P ( 


-L -L 

f ^ u x +a) 9 yn z dz + J" ^ 2 “x _a) 

0 0 


s z dz 
xn 


- r <■ 

0 

5 * 


( u x + “ y) 9 n 2 dz ^ 


sin 0 L L cos 0 L 


since \ zsinB z dz » - — r- 

6 n 


0 


f 3 < 6 n> 


-L 

f zcos0 z dz 
n 


I 


zsinh0 z dz « - 

n 0 


cos 0 L L sin 8 L 

— r — + — s — 2 !--*.<•> 

S 2 °2 6 2 4 n 

n n 


L cosh 0_ L sinh 0 L 


— + ^ - f 5 (8 n> 


0 


n 


and 


I 


zcosh 8 z dz 
n 


L sinh 6 L cosh 6 L , 

- 2 s — + - f (6 ) 

e 2 e 2 « n 

n n 


8 


sin 0 L L cos 0 L 

* n * p K' "pT + 0 2 ) [ (u 2 <») x -d) y ) A ln +(ajzwy+ux) A 2n ] 


L sin 0 L cos 0 L-l 


+ T 1 i y >B la +( “z“y^x )B 2n 1 
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slnh 6 L L cosh 6 L 

+ ( -y-S g — : )t< “z" I[ -“), )C l „ +< “ z “y“x )C 2n 1 

P_ n 

n 

L sinh 6 L cosh 6 L-l 

+ (t - n — » n ) I (m u» -a) )D, +(oj m +d) )D. ] 

6 a 2 a x — y In z y x In 

n p — 


1 1 

L cos 0 L sin 0 L 0 9 

-( ^ 


8 6 
n., n, 

L sin 0 L cos 0 L-l 


2 . 2 , 


-< ><•; + -; )b 3« 


0 . 


6 


n 


or 


M 

* r <£ 3 ( 8 0 ) K “* “* - “V A ln + ( Vr + “x 5 4 2 a ] 


+ W [ “z“x -“ 7 )B ln + (,,; z“ T + “x )B 2n ] + f 5 (S n ) [ “z“x ' “P C ln 


+ z“y + % )C 2n ] * V 6)t “z“x * “y> D l„ + <“*“ y * “ x> D 2n ] 
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oos <pooa 9 


sirup 


sLnd cos<f> 


sinpsin^ coefsin <PcoeQ cosiJ>cos<fr cos^ain^ai n 8 + si n ft cos 8 
£rfntsdni|>cos 0 + oos^sinQ simp cos<P ooa ij/cosQ-sin^ aln <ft sin 9 



COS COS 0 

sin <p 

-sine cos (p 


sin ^ sin 9- cos^sin^ooee siity sin cose + cosi|/sin<j> 
eoBipcoa <p ■ - sln<jicos<j> 

sin ^ cos 9+ oostj; sin 41 sin e • oos^cosS- sin^ sin4> sin 9 


Assuming here the intrinsic frame^^ of reference to 


coincide with 

the 

orbit frame, 1. 

e . 

X » 

0. 



M u 

M 12 

M 13 


0 


* M 13 ' 

and M°(? o ).„2 

M 21 

M 22 

M 23 


0 - 

2 

% 

zM 23 


M 31 

. 

M 32 

M 33 


z 


N 

S 

CJ 

u 

1 — 


where M 13 * 3u) o (S6CdGpGl>- S 2 9S<J>SiJiC $ 

^23 ■ 3 u)q [SeCeS<f(l- 2 C 2 i|)) + S t|/ C i|< (S 2 $S 2 9 - C 2 0)] 
M33 - m 2 {-3(s^sesifi - cec ip ) 2 + 1} 


and 
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Therefore 

*o f* • M<r 0 ) <■“ ■» j < s xo* M 13 + V ”23 + z9 n *33 1 iz 
M 0 

*n • r“ (f 3 ® n> (i ln «13 ’ A 2n «23> + W (B ln M 13 + B 2n «23> 

+f 5 ( 8 n> «1„ M 13 + °2„ «23> + W <“ln M 13 + D 2n M 23 ) 

+ W A 3n Hj 3 + f 4 (s n ) B 3n Mjj 

the generic mode equations can be recast as: 

A n + “n A n + £ ' ' “y * V 1 f 3 (s n )A to + f * (e n )B ln + Win 

* W “to 1 + ( “z "j * % - «23> 1*3® n > A 2 B * W“ta + *5 (V<2. 

+ ' W “2. 1 + -“x -“J> 1£ 3< 6 ’n> ^ + V 8 ’o> V * 

• • 

P W-° W 0 + + 3 yn *f-» 
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The scope of this study is limited to the stability 
analysis of the SCOLE system and the derivation of a 
control strategy for its motion about the nominal 
equilibrium position as derived in Chapter IV. The 
gravity-gradient forces acting on the system will be also 
calculated in that configuration. 

Thus, in what follows, the Euler angles ip , 0, and<j> 
will be replaced by rij » r^* and with tj; ■ ^ e q + nj , 9 - 
9 eq + n 2* ^eq + n 3 * After linearizing the different 
terms appearing in the generic mode equations, they can be 
rewritten, for each of the four modes included in this 
study, as: 

*n ♦ % A „ + T ( - ( S 2 + 3 “0 "2> t W A ln + W <> ln 
+ f 5 ^ Sn + f 6 “in 1 + (n l' 2 V? 4 “ *2n + £ 4 \ 

* f 5 (f, -J S. ■ * f 6 V V + \ ' 2 I h V S, + f 4 V 

■ F x WO * Vfc. <-«W + VV H), 3 + V“ > v 2 + 3 jn 


s 

/ 
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In Chapter V, these equations (n » 1, 2, 3, and 4) 
will be added to the equations describing the rotational 
motion of the SCOLE to obtain a mathematical model of the 
SCOLE orbiting configuration. A modified version of 
equations III. 22 - III. 24 will then be used. 



